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SUMMARY 


In this report, the work carried out at the University of Kansas 
Flight Research Laboratory (KU-FRL) to measure the transmission loss 
characteristics of panels using the acoustic intensity technique is 
preriented. The report describes the theoretical formulation, installation 
of hardware, modifications to the test facility, and development of 
computer programs and test procedures. A listing of all the programs 
is also provided. 

The initial test results indicate that the acoustic intensity 
technique can be easily adapted at the KU-FRL test facility to measure 
transmission loss characteristics of panels. Use of this method will 
give average transmission loss values. The fixtures developed to 
position the microphones along the grid points are very useful in 
plotting the intensity maps of vibrating panels. 

Based on the experience gained so far, several improvements to 
the test facility and test procedures are also identified. 
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CHAPTER 1 




INTRODUCTION 

This report iw a continuation of the documentation of the research 
accomplished under continuing NASA Cooperative Agreement NCCI-6. The 
progress of the research accomplished during the period November 1982 
through April 1983 of the previous project year (May 1, 1982, through 
April 30, 1983) was included in the previous report, KU-FRL-417-21 
(Reference 1) . 

The present report covers the period from May 1, 1983, through 
November 30, 1983. In this reporting period, the measurement of 
transmission loss characteristics of panels using acoustic intensity 
techniques is being investigated. The characteristics and the 
limitations of the present measurement techniques are described in 
Reference 2. Most of the limitations mentioned in Reference 2 are 
due to the small size of source and receiver sections of this test 
facility and the use of acoustic pressure levels instead of acoustic 
sound power levels as a measure of sound power. Also, in the present 
method, the sound pressure levels are measured at only one location. 
Even though this location had been chosen after a careful experimental 
study, it is possible that this location may not be ideal for some 
cases. Measurement of the sound power by the integration of the 
acoustic intensity levels over the entire panel will eliminate a 
few of these limitations. The direct measurement of the .acoustic 
intensity has now been made possible by the development of the two- 
microphone, cross-spectral method. This report describes the adoption 
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of chls measurement technique at this test facility to measure panel 
transmission loss values. 

Theoretical developments for the calculation of the acoustic 
intensity from the pressure measurements by two microphones separated 
by a known distance is given in Chapter In the same chapter, some 
of the limitations of this method and ways to reduce some of the errors 
encountered are also described. The present test set-up had to be 
changed to introduce the intensity method at this test facility. The 
modified teat set-up is presented in Chapter 3. Also presented in this 
chapter is the description of the computer programs and the modified 
test procedures. A typical test result from this test facility obtained 
using the acoustic intensity technique is given in Chapter 4. The 
comparison of the present test results with the previous results, and 
the modifications that will be required to Improve test time and 
memory requirements are presented in ’Chapter 5 . 
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CHAPTER 2 


THEORETICAL ANALYSIS 


2.1 ACOUSTIC INTENSITY 


The acoustic Intensity at any point is defined as the rate of 
acoustic energy flow across a surface of unit area (Reference 3). 
By definition: 


I 


r«inst 



( 2 . 1 ). 


This energy flux, 5^ , Is equal to the amount of work done upon 

the area A in the direction r due to the total force, F^; i.e., 

• 6r » p^5^ • <5r (2.2), 

r 

where p. is the total pressure comprising the ambient pressure p^ 

C a 

and the sound (perturbed) pressure This gives 


r.inst 


p^u + pu 


(2.3), 


where u ■ dr/dt is the particle velocity in the direction r. Both 
the sound pressure and the particle velocity are functions of spatial 
coordinates and time. For sinusoidal processes, the time-averaged 
value of the first term is zero if the averaging time is an integral 
number of half periods. For other processes, it well be zero if the 
averaging time is sufficiently long. If the processes are stationary 
randoms same result can be obtained by 

E{I } - E{p^ • u^} + E{pu^} 

17 3 i t 


= p^E{u^} + E{pu^} 


= p_u + E{pu } (2.4) . 

^mean 
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If the mean flow Is zero, then 

E{lj.} - E{pu^} (2.5). 

Direct measurement of intensity using pressure-velocity product 
has proved very difficult in field conditions (Reference 4). An 
Indirect measurement, wherein two microphones are used to measure the 
acoustic intensity, has gained wide attention in recent years 
(Reference 4). In the next section, equations required for the 
measurement of acoustic intensity using this method will be derived. 


2.2 ESTIMATION OF ACOUSTIC INTENSITY USING TWO-MICROPHONE METHOD 

With zero mean flow of the medium, the time-averaged intensity is 
given by Equation (2.5), For ease of calculation, let us consider both 
p(r,t) and u(r,t) to be stationary random processes. Fourier transforms 
of nttciilonary random processes exist if their autocorrelations and cross 
correlation are aperiodic (Reference 6), In such cases the Fourier 
transforms of p(r,t) and u(r,t) are defined as 

00 

P(r,to) » - / p(r,t)e"^**^^dt (2.6) 

«p09 

00 

U(r,u) » - / u(r,t)e^**^^dt (2.7). 

From Euler’s equation (Reference 3), the relationship between 
the particle acceleration and the pressure is obtained as 

3 

P ■ -grad p (2.8). 
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In one direction, namely r, 


3t 


_i£ 

9r 


(?..9). 


In subsequent discussions, it is assumed that the particle velocity 
is in the direction r, and hence the subscript r will be dropped. 
The particle velocity is obtained by integrating Equation (2,9): 


u - - ~ dt (2.10). 

p ^ 9r 

To measure intensity using two microphones, an Intensity measurement 

apparatus as shown in Figure 2.1 is used. In practice, the pressure at 

the center of closely spaced points A and B can be approximated by 

taking the mean of p^ and pg. The pressure gradient, to a first order, 

can be calculated by dividing the difference in pressures at p. and p„ 

A 15 

by the separation distance, dr. These approximations give the following 
estimates for p(r,t) and u(r,t): 


p(r,t) = j (P^(r,t) 




1 *• 

O 

This approximation can be considered valid as long as the 
separation is small compared to the wavelength, X (Reference 5). 
Following Laplace Transform procedures, the time integral of the 
transform can be replaced by 


( 2 . 11 ). 


^[/udt] = (2.12). 

Reference 5 states that even though this procedure is mathematically 
incorrect, it gives valid results in practice. Hence Fourier transforms 
of p(r,t) and u(r,t) can be written as 
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OF POOR QUALITY, 


Dlrecclon of 
Acoustic Intensity 



Microphone B 


Microphone A 


Figure 2.1: Acoustic Intensity Measurement Apparatus 
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P(r,w) « j {P^(r,o)) + Pg(r,w)} 


U(r,w) 


From Equation (2,5): 


jw6rp B 


{P_(r,u) - P.(r,w)} 


^r,av “ E{p(r,t)u(r,t)} 


(2.13) 

(2.14) . 

(2.15) . 


Both p and u are functions of the spatial coordinates of r. The 
cross correlation function of p and u is defined by (Reference 6); 

R (tj^t2) “ E{p(tj,)u(t2) } (2.16), 


Because of the stationarity * this equation can be written as 

Rpu('T) * E{p(t)u(t + t)} 


(2.17). 


At T * 0, 


Rpu(O) = E{p(t)u(t)} 


(2.18). 


The right hand side of the equation is equal to the averaged intensity. 
Therefore, 


I ^ R (0) 
r,av pu' 


(2.19). 


By definition, the cross-spectrum of these two processes is given by 
(Reference 6) 

pa 

Sp^(r,to) = / Rp^(T)e”^^^dT (2.20), 


and its inverse Fourier transform is 


R (t) = i / S (o))e^‘^\iLJ 
pu ' ' 2tt pu ' 


( 2 . 21 ), 


With T = 0, 




( 2 . 22 ) 
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and 

R (0) - R (0). 
pu ' up' 

If the Fourier transform of p(t) and u(t) exist, the cross spectrum 

can be written as (Reference 6) 

S - E{P(f)U*(f)} (2.23). 

^u 

Substituting the values for P(f) and U(f) from Equations (2.13) and 
(2.14), 

E{PU*} * E{[i (P^ + P^)lt- 3^ (Pg - P^))*) (2.24)- 

Simplifying this equation, 

^'W> - ®<Va*H <2.25). 

By definition, 

E{PgPg*} = Power spectrum of pressure at ^ “ ^AA» 

E{P^P^*) “ Power spectrum of pressure at A = Sgg, 

E{P^Pb*} * Cross power spectrum between pressure at B and A * S^, 

E^^B^A*^ =• Cross power spectrum between pressure at B and A * 

Because Sg^ - S^*, 

^ ^^AB ~ ^BA^ “ +2Im(Sg^). 

Substituting these relations in the equation, 

= 2 ^ <®B - =AA> 2Im(S3^)} (2.26). 

If the cross correlation is real, which normally is the case, 
the real part of the cross spectrum wall be even and the imaginary 
part of the cross spectrum will be odd. Hence, when integrated from 
- 00 to + “, the odd part integrates to zero. Using only the real part. 
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I - 

r.av 

00 




(2.27) 


Fourier analyzers use only one-sided spectrum. The values on 
the positive frequency side are doubled to keep the energy the same. 
One-sided cross spectrum is normally denoted by 


X = -/ — j — Im(G_.)dto 
r,av 0 BA' 


(2.28) 


The negative sign in the equation can be avoided if the microphone 
closest to the source is connected to channel B of the analyzer (see 
Equation 2.11). For this case, the intensity can be written as 


I = / Im(G.^)dw 

r,av - (jj5rp AB 


(2.29) 


In practice, the digital form of the estimate will be used: 

, N/2 IraG.^(nAf) 

I » ^ 

r.av p6r ", nAf 
n=l 


(2.30) 


where Af is the calculatioi’i bandwidth and N is the block size of the 
analyzer. Intensity as a function of frequency is 


I^(nAf) 


^ Im G^g(nAf) 
p5r nAf 


(2.31) 


2.3 LIMITATIONS 


References 3 through 5 discuss the inherent limitations of the 
two-microphone cross-spectral method to estimate the acoustic intensity 
The limitations arise due to two types of error that occur: 
a) a systematic error and b) a statistical error. The systematic error 


is due to the finite difference approximation used in the formulation 



of acoustic intensity. The statistical errors are due to the random 
source excitation and other random variations in measurement. In 
addition there are some more limitations that are specific to the 
KU-FRL acoustic test facility. All these limitations are discussed 
below. 


2.3.1 High Frequency Limitation 


At the KU-FRL acoustic test facility there are two possible sources 
of error in the high frequency region. The first limitation is due to 
the finite difference approximation for pressure and pressure gradient. 
This produces a systematic error in the estimation of these two 
quantities. The approximations used are (Equation 2.11) 


Pa + Pb 


(2.32) 


3 £ 

3r 


Pb - Pa 

iSr 


(2.33). 


By the mean value theorem, these approximations tend to the actual 
values only when the separation distance tends to zero. Otherwise, 
they produce a systematic error in the entire frequency range. However, 
the error is most severe in the high frequency range. For a plane 
sinusoidal wave, the estimate of the intensity, using this approximate 
method, is related to the actual intensity by (Reference 5) 

I _ sin(k6r) 

I ” k6r 


where = actual intensity, 

“ calculated intensity, 


k ■ wave number (w/ c ) . 
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(Sin x/x) tends to 1 when x tends to zero. Otherwise, it is less 
than 1. Hence at high frequency (high h) and large separation distance, 
the acoustic intensity id.ll be underestimated. At the KU-FRL acoustic 
test facility, this is minimized by limiting the separation distance 
to 25 mm (1'*) at frequencies above 500 hz. 

The second limitation is due to the band pass characteristics of 
microphones. Because low frequency noise reduction is the major concern 
in aircraft noise reduction, microphones with higher sensitivity in this 
region are preferred. Because only the low frequency region is important 
in aircraft noise control applications, 1/2” B&K microphones were chosen 
for the measurement of transmission loss of panels at the KU-FRL acoustic 
test facility. These microphones are accurate only up to 4000 hz. With 
1” separation and up to 400 hz, the error due to the approximation will 
be less than 3 dB for a plane with sinusoidal wave. However, because 
this is a systematic error, similar error occurs with and without the 
panel. Hence, when the transmission loss is calculated, these errors 
tend to cancel each other out. 

2.3.2 Low Frequency Limitations 

According to Reference 4, there is no evidence of any low frequency 
limit due to the approximation errors. Reference 5 shows that the 
estimation of the particle velocity results in the estimation of the 
phase angle difference between the two microphones. The term '’k,5r” in 
Equation (2,34) is the phase difference between the microphones. This 
term is very small at low frequencies because k is small. Hence, at 
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low frequencies, the measurement error of the phase angle becomes 
significant. The measurement error is due to the channel mismatch 
between the two microphone channels. This error can be eliminated 
(or reduced) either by using phase-matched microphones or by correcting 
for the difference in the phase angles when both the microphone channels 
are exposed to the same sound field. While the use of phase-matched 
microphones will make measurement easier, it cannot account for the 
phase mismatch in the rest of the measurement channels (like signal 
amplifier, etc,). Because of this, a phase calibration procedure 
is being adopted at the KU-FRL acoustic test facility. These procedures 
are discussed in detail in Chapter 3. 


2*3.3 Near Field Limitations 


The third limitation occurs when this method is used in cases 
where the intensity changes rapidly along the probe. When this occurs, 
the intensity is very different at the two microphone locations. A 
similar situation arises when the measurements are made in near field. 
Several expressions have been derived to estimate the effect of 
near field for simple sources such as monopole, dipole, and quadrupole. 
Table 2.1, taken from Reference 5, gives the following criteria for 
limiting this error. 


Source type 
Monopole 
Dipole 
Quadrupole 


Proximity error less than 
1 dB if source is away by 

>1.1 6r 

> 1.6 5r 

> 2.3 6r. 
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While these results will not be valid for a complex source such as 
a thin panel, they do provide some guidance in using the acoustic 
intensity techniques near the sound sources. 

2.3.4 Limitations Due to Statistical Errors 

Because of the random excitation, an estimate of (f ) is made. 

This estimation gives an additional error due to the variance of the , 

quantity being measured. Reference 7 gives the normalized random 

1/2 

error, e(I) ■ (Var(I) /I), in this type of measurement as 

e(I) - (n^)"^''^[l/Y^ + cot2^,^(l - (2.35), 

where n^ is the number of ensemble averages for cross spectrum, and 
Y^ is the coherence between the acoustic pressure at the two measurement 
points. 

As can be seen, the statistical error can be minimized by selecting 
a large number of ensemble averages and making sure that the coherence 
level is high. Since the tests are conducted inside a closed cavity 
where no other sources exist, the measured coherence values are 
normally vary high. In the KU~FRL acoustic test facility an ensemble 
average of 256 and acceptable coherence of above .8 are used. For an 
assumed phase difference of .18 rad, with these values for ensemble 
averages and coherence, the statistical error (e(I)) will be less than 
.194. For a plane wave, a phase angle difference of .18 rad corresponds 
to 100 Hz at 4" microphone separation. 


2.4 CORRECTIONS FOR PHASE MISMATCH 


As discussed in Section 2,3, phase mismatch between the two 
microphones can be minimized either by using phase-matched microphones 
or by correcting for the error. One of the disadvantages of using the 
phase-matched microphone is the error due to phase mismatch of the 
rest of the measurement channel cannot be corrected. At times these 
errors may become si^if Leant. Hence at the KU-FRL acoustic test 
facility, phase correction by prior calibration of microphones is 
used. A literature search was conducted. Based on the results, four 
promising methods were chosen (References 4, 5, and 8). 


2,4,1 Phase Angle Correction 


In this method the phase difference between the two measurement 
channels (including microphones) is measured when the microphones are 
subjected to the same sound field. The phase angles of the cross 
spectrum measured during the intensity tests are corrected for this 
difference. The magnitude correction is done separately. If the 
same sound field is applied to both the microphones, shown in Figure 
2.2, the measured cross spectrum is given by 

S._ - • H.* • H^ (2.36), 

AB Pj^P2 A 3 


where S is the cross spectrum of the sound field at the position 

P1P2 

of the two microphones, is the measured cross spectrum, and 
and Hq are the transfer functions of the two measuring channels. 

The phase angle of the measurement channels is the phase angle of the 
transfer function. 
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This is one of the methods chosen at the KU-FRL to correct for 
the phase angle difference. This method is useful at low frequencies. 
The exact realization is discussed in the next chapter. The magnitude 
calibration is done separately using B&K Plstonphone 4220. 


2.4.2 Transfer Function Method 

Reference 5 shows that when two microphones are exposed to the 
sound field, both magnitude and phase correction for channel mismatch 
can be done using the relation: 


PiP 


'AB 




(2.37), 


where H._ is the transfer function between the measurement channels. 
AB 

Since this method is very similar to the previous method, this was 
not tried. 


2,3.4 Microphone Switching Method 

Chung, et al. (Reference 4), originally proposed this method for 
correcting phase mismatch. In this method, tests are done twice. 

Tests are first performed with the microphones in normal locations; 
tests are then repeated with the microphones interchanged. Under 
these conditions Reference 4 gives the actual cross spectrum as 

Im = UG^ • G^]^^^}/p6rw • 1 h^| • iH^j (2.38), 
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where * cross spectrum between microphones, 

S 

■ cross spectrum with microphones switched, 

|h^|, |Hg[ * gain factors, microphones A and 
In this method every test has to be done twice; also, therefore, the 
test section has to be opened for every measurement. For these reasons 
this method is not being used at the KU-FRL acoustic test facility. 


2.4.4 Modified Microphone Switching Technique 


This method is a combination of the transfer function method and 
the microphone switching method. In this method, before the start of 
the tests, the microphones are exposed to a sound field and the* cross 
spectrum is measured. Now the microphones are switched, the 

measuring system is exposed to the same sound field, and once again 

c 

the cross spectrum is measured (G^)* From Reference 8, we get 


M 


4 . 




(2.39), 


where (|i is the phas-e angle between the measurement channel. By 
assuming that the magnitudes are the same, the complex root computation 
is avoided. The phase angle is calculated by dividing the phase angle 
of the cross spectral division by 2. 

This error is used to correct the measured intensity values 
during the actual tests. The implementation of this method at the 
KU-FRL acoustic test facility is discussed in the next chapter. 

The advantages of this method are 1) the microphones need not be 
exposed to the same sound field, 2) tests need not be performed twice. 


3) the method is valid even at high frequencies. The only requirement 
is that the sound field should be stationary. 


CHAPTER 3 


EXFERII’IENTAL SET-UP 

* 

3.1 HARDWARE DESCRIPTION 

3.1.1 General Test Set-Up Description 

The general arrangement of the acoustic Intensity test set-up 

Is shown schematically In Figure 3.1. The system shown was designed 

to take and process data as quickly and efficiently as possible. Since 

each TL test requires 324 Intensity spectra at 402 frequency values 

each (325 spectra ■ 81 points for high and low frequency tests for both 

the source and receiver side) , the need for speed In data processing 

and efficiency In data storage 1 *comes obvious. The operation of the 

system is described below. 

* 

The heart of the system is the Nicolet 660B dual channel FFT 
analyzer. The analyzer provides temporary data storage and performs 
all required FFT calculations. It is controlled by a Zenith Z-100 
microcomputer which provides data reduction and permanent data storage 
capability. The 660B and Z-100 ar- linked through their respective 
RS-232C ports at a 9,600 baud rate. At present, the communication 
software used to transfer data from the 660B to the Z-100 is written 
in a high-level language and represents one of the limitations to the 
speed of testing. The development of an assembly language communication 
program (a future project) will, however, decrease the test time. 

In addition to its data acquisition role, the Nicolet 660B also 
provides the excitation signal that drives the speakers in the Beranek 
tube. This excitation signal is a band-limited binary white hpise 
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output from the analyzer’s rear panel. It is passed through a TAPCO 
2200 equalizer for purposed of modifying the speaker Inputs to achieve 
a flat speaker output. The equalizer output Is gained up through a 
Crovn D-150 power amplifier to drive the nine Altec 405-8G loudspeakers. 
It may, however, be necessary to Insert a high pass filter between the 
analyzer and Lhe equalizer when testing panels with large transmission 
losses. This would be required to avoid overloading the analyzer 
inputs in the low frequency range when attempting to gain up the 
microphone outputs in the high frequency range. As yet, this had not 
become necessary. 

Two B&K 4165 microphones with B&K 2618 preamps are positioned in 
the Beranek tube by the microphone positioning device (MPD — described 
in the next section) . The microphone preamplifier outputs are fed 
into the two channels of the 660B FFT analyzer (although tests involving 
panels with very high transmission losses may require additional 
amplification of microphone signals — such as the KU-FRL Nagra SJS tape 
recorder— between the microphone power supply and the analyzer) . From 
the analyzer, the cross spectrum of the two microphones is transferred 
to the Z-100 microcomputer where it is stored on 5 1/4 inch disks. 

Data transferred to the Z-100 are cataloged in files by microphone 
location, analysis (frequency) range, and source or receiver spectra 
so that batch processing of data is simplified. Data reduction routines 
are run on the Z-100 to generate point intensity values and overall 
panel transmission loss. These data are transferred to a Digital MING 
computer for plotting on a Hewlett-Packard Model 7225B X-Y plotter. 
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3.1.2 Microphone Positioning Device (MPD) Description 


The microphone positioning device was designed and built at the 
KU-FRL for the purpose of accurately positioning the microphones within 
the Beranek tube. The design requirements specified that the MPD be 
able to position two microphones anywhere in a 16 inch by 16 inch plane 
parallel to and directly behind the test panel without opening the 
tube. Movement of the mikes had to be done easily and accurately from 
the outside. In addition, provisions for varying the spacing between 
the microphones had to be made, and "blockage" due to the device 
(interference with the sound paths within the tube) had to be kept to 
a minimum. A later requirement that the MPD provide for easy time-area 
averaging could not be met, as construction had already progressed too 
far at that time, 

* 

The MPD is shown in Figure 3,2. It is an extension tube constructed 
of particle board into which the positioning mechanism is built. 

Vertical and horizontal motion is provided by a system of cross beams. 

A Lucite block is attached to the vertical and horizontal beams at their 
intersection and is allowed to slide freely on both. The block is 
therefore constrained by the cross beams (guide rods) such that when 
the rods are moved, the Lucite block maintains its porition at their 
intersection. The microphones are attached to the Lucite block through 
an aluminum beam protruding from it (see Detail A). The microphones 
can be positioned at different locations along the beam to provide for 
different microphone spacings. 
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Figure 3,2: Description of Microphone Positioning Device 






Detail A 




I'he guide rods in the MPD are controlled externally by a cable 
and manual crank system. Position information is displayed on scales 
by a secondary cable system driven off the cranks. 

The MPD operates smoothly and positions the microphones with 
reasonable accuracy. However, due to interference of the microphone 
cables with the bottom of the MPD at low positions, it is not possible 
to cover the entire 16 inch by 16 inch sweep area. The solution to 
this problem is to turn the microphones face down when they are 
positioned near the bottom of the MPD, This, however, requires that 
the Efisranek tube be opened midway through a test. While this is not 
a significant problem, it Increased testing time. 

3.2 SOFTWARE DEVELOPMENT 

Because of the large amount of data that will have to be processed, 
using this method, the computer program had to be split into many 
subparts before it could be handled by the Z-100 computer. Depending 
upon the ease of programming and the amount of calculations involved, 
either Fortran or Basic language was chosen to write these programs. • 

The Flow diagram shown in Figure 3.3 describes the steps Involved. The 
individual steps and the relevant equations are described in subsequent 
sections. 

3.2.1 Magnitude Calibration 

A B&K "Pistonphone” is used to calibrate the microphones. Because 
the 660B outputs unsealed values, the actual output from calibration 



Figure 3.3: Test and Analysis Flow Diagram 
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tests Is a function not only of the pressure but also of the input max 
amplitude setting and number of ensemble averages. In converting the 
output of the 660B to the actual BNC input volt level and then to 
pressure, these. two additional variables will have to be considered. 

The B&K 4220 Plstonphone outputs calibrated sound pressure level 124 dB 
re 20 micro pascals at 250 hz. Hence, 


20 log(-£2i) » 124 dB 
^ref 


20 ^ 

10 p. 


(3.1), 


'cal ‘^ref 

where P^^i “ pressure corresponding to pressure level of 124 dfi 
Pj.g£ =• reference pressure (20 micro pascals) . 

At a given input channe], maximum amplitude setting and for a given 
number of ensemble average, at any cell location of 660B, the pressure 
(p^) will be proportional to the value output by the 660B (v^) . 

Pi » 


or 


Kv^ 


(3.2), 


wher@ K. is the calibration constant. The Plstonphone outputs 125 dB 
sound level at 250 hz. There is a small tolerance about 250 hz. 

Also, spectral leakage always exists in digital signal processing. 
Whenever the energy is concentrated at a discrete frequency which is 
in between two adjacent cell (filter) locations, the energy is smeared 
across the neighboring cells. See Reference 9 for discussion on 
spectral leakage. In order to minimize the effect of spectral leakage 
during calibration, the power-spectral values of three adjacent cells 
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on either side are summed to obtain the total energy, the calibrated 


pressure can be equated to 


lg+3 


E , Pi" • kEvi^ 


l-ig-3 


cal 


io+3 


K E V 
i-io-3 


(3.3) 


(3.4), 


where i^ is the filter location corresponding to 250 h^, and is 
the value output by the 660B at a given maximum amplitude setting 
and for a given number of ensemble averages, the calibration constant 
K can be calculated. This needs to be done for both channels. The 
functional relationship between the output and the ensemble averages 
and the maximum amplitude Setting is given in Reference 10. Based 
on these relationships, the relationship between the true value and 
the value output from the analyzer 660B during any one test was 
derived as follows. 

RMS spectrum of channel A: 


TV - K. • V. • 7^ • /nTn. 
A t A c t 

c 

Power spectrum of channel A: 


(3.5); 


A N 

TV • K • V • )2 • — £• 

’'A ^ N. 
c t 


(3.6); 


Cross spectrum: 


'^Vt • h.t'> 


N 


• \,c> ‘'c 


(3.7); 
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where TV ■ true value, 

V ■ value output, 

A • maximum Input amplitude setting, 

N ■ number of ensemble averages, 

K « calibration constants defined in Equation (3.3), 
and the subscripts t, c, A and B correspond to test, calibration, 
channel A and channel B, respectively. These relationships were 
confirmed by experimentation. They are used in obtaining calibration 
constants. The actual test and analysis procedure developed, based 
on the above equations, is described in Appendix A. The listings of 
programs PSP660 and MAGCAL, used for the determination of magnitude 
calibration constants, are given in Appendix B. The output from these 
programs are stored in a file named CALDAT.DAT. It stores calibration 
factor, number of averages, and maximum amplitude setting for both 
channels. This file is accessed by other routines to convert the test 
values into true values. 

3.2.2 Phase Calibration 

As described in Chapter 2, two different calibration techniques 
are used at the KU-FRL acoustic test facility. Method 1 calculates 
the phase angle difference between the two microphone channels when 
both the microphones are exposed to the same field. Method 2 uses the 
modified transfer function method described in Chapter 2. 
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3. 2. 2.1 Method 1 


In this method both the microphones are exposed to the same field, 
and any difference in the phase angle me -"tred is due to the difference 
in the channels. Subsequent tests can then be corrected for this 
difference in phase angle. Figure 3.4 shows the schematic diagram for 
the microphone phase calibration system. In this method, the two 
microphones are inserted into a long tube with faces of the microphone 
parallel. A random noise is generated at the other end. Hence both 
the microphones are exposed to the same sound field. Only the cavity 
resonance effects affect the actual sound field incident at the 
microphone. By the tube diameter of two inches, the fundamental 
circumferential resonance frequency is made to occur at a frequency 
greater than 3000 hz, which is the maximum frequency of interest. 

Thus the effect of circumferential resonance frequency is avoided. 

The effect of longitudinal resonance frequency could not be eliminated 
fully, but it is minimized by having absorptive fiberglass materials 
on the ends of the calibration tube. 

During the initial determination of the phase angles, it was 
noticed that a certain amount of scatter was unavoidable in the phase 
angle differences measured. Since this scatter may affect the results 
during daily calibration, a statistical approach was taken to minimize 
the effects of this scatter. It was decided to perform tests many 
times to cover the entire range of parameters that cannot be controlled 
exactly during any test. These parameters involve the humidity, 
temperature, amount of time the calibration speaker has been on, etc. 
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Thirty tests were conducted to cover the range of variables. 

A mean of the results of these thirty tests can be considered to be 
a good estimate of the mean of the population of all possible phase 
angle measurements (see Reference 11), However, thirty calibration 
tests every day to cover all possible random combinations is not 
practicable. Hence it was decided to use significance testing 
to obtain acceptable calibration values. In this procedure, 
the population mean and standard deviation are first determined only 
once. Thereafter, only a small number of tests need to be done every 
day. The mean values of these tests are compared with the population 
mean values, and the significance tests are used to accept or reject 
the new values. 

An estimate of population mean can be obtained by taking a mean 
of a large number of tests. If t'K. ■•umber of samples is greater than 
thirty, it can be assumed that the mean and the standard deviation of 
the sample are equal to the mean and the standard deviation of the 
population (Reference 11). Hence, thirty tests that are conducted 
in the beginning of a test series can be assumed to be a very good 
estimate of the population mean and the standard deviation. Daily 
calibration values are then compared with these values for acceptability. 

Acceptability or rejection of a test sample (in this case daily 
calibration results) is based on the significance testing. The aim 
is to minimize both type I and type II errors, A type I error is 
committed if the null hypothesis is rejected when it is true, A 
type II error is committed if the null hypothesis is accepted when 
it is not true. In this case, the null hypothesis will be that the 
sample mean (x) is equal to the population mean (Pq) » l.e., 
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(both means are the same) 


(3.8) 


HO: X ■ Uq 

The alternate hypothesis is given by 

HI: X (the mean of the sample is not equal 

to the population mean) (3*9). 

The test statistic is given by 



(3.10), 


where TI 2 Is the number of averages of the test sample (during dally 
calibration) and 0 is the standard deviation. 

HO is true HO is false 

Accept HO Correct decision Type II error 

Reject HO Type I error Correct decision 

Specifically in this case, while committing type I error can be 
tolerated, committing type II error should be avoided. The probability 
of committing type II is denoted by The probability of committing 

type I error is denoted by a. This is also known as the level of 
significance. When the alternate hypothesis is nonspecific, as in this 
case, it is not possible to compute the probability of computing type II 
error (Reference 11). However, with higher sample size, both a and B 
can be reduced. Reference 11 also gives the following equation for 
the two-tailed test to obtain the power (1-3) for a specified 
alternative as 


"2 " 




(3.11), 


where 02 is the number of observations required, a is the standard 
deviation, and 6 is the difference between the sample mean and the 
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is 1.96 for 


population mean. For .05 level of significance (a), 
normal distribution, and for .05 probability of computing type II 
error (0), z •• 1.645. Using Equation (3.11) as a guide and by trial 
and error, n » 5 was observed to be adequate for our calculations. 

These equations have been modelled into the computer program. 

At the beginning of a series of tests, the calibration is performed 
30 times, varying the uncontrollable parameters (such as temperature, 
humidity, etc.) as much as possible. Then a basic program, STAT.BAS, 
is run. These tests are performed once for low frequency range and 
again for high frequency range. The outputs (the population mean and 
the confidence interval at 95% confidence level) are stored into two 
files. For t distribution, the standard deviation (a) and confidence 
interval (c) are related at 95% confidence level by 

c * /tvT 

^ “ 2.02 (3.11), 

where n^^ is the number of tests. 

During the day of the tests, calibration is done only five times. 
The analysis program, CALII.EAE, is run to perform the significance 
tests. During the running of this test, the operator xd.ll be asked 
whether to accept or reject the calibration if the z statistic 
(Equation 3.10) exceeds 1,96. At this z statistic, the level of 
significance is 5% and the power is about 95%. Very high values of 
z statistic should be rejected and the calibration redone. This has 
to be done for both frequency ranges. The output file from this 
program is called CALII.DLO, or CALII.DHI. These files contain the 
cosine and sine of the phase angle correction at each filler. These 
files are accessed by other routines to correct for phase angles. 
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3. 2, 2. 2 Method 2 


The second method for phase correction uses the modified microphone 
switching technique described in Section 2.4.4. Equation (2.39) is used 
to obtain the correction. In this method, the tests are done only at 
the beginning of a test. First, the microphones are clamped in normal 
location in the MPD and the cross spectrum is measured. Then the 
microphones are switched and the switched cross spectrum is measured. 

From these two cross spectra, the phase angle correction as a function 
of frequency is obtained using Equation 2.39. The test procedure is 
described in Appendix A; and the listing of the program INTCAL, which 
performs the calculations, is' given in Appendix B. 

3.2.3 Intensity Tests 

The test procedure for measuring acoustic intensity values at 
the KU-FRL acoustic test facility is given in Appendix A. The intensity 
is calculated from the measured cross-spectral values by Equation (2.31). 
The program INTSTY performs this calculation. It also performs relevant 
magnitude and phase corrections. At present, the intensity values are 
calculated at 81 grid points on an 18-inch-by-18-inch cross sectional 
area. These intensity values can be made use of to plot either an 
intensity map or transmission loss. The relevant programs are identified 
in Figure 3.3. The listing of programs is given in Appendix B. 


3.2.4 Plotting 


r>*=» 


' I 

%i 


il : 


At present, no plotter is available for the Z-100 microcomputer. 
Hence the plotting is done- on a MINC-11 computer with a HP 7225B 
plotter, 'fherefore the data values have to be transferred from the 
Z-100 to the MINC-11 computer. This is achieved through modem 
communication. The data values are then plotted using plot programs 
ZPLOT and PLOT (for transmission loss and intensity maps, respectively). 
The listings of these plot programs are given in Appendix B. 
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CHAPTER 4 


TEST RESULTS 

This chapter describes the tests conducted to check out the acoustic 
intensity procedures developed at the KU-FRL acoustic test facility. At 
the time of preparation of this progress report » this activity has not 
yet been completed. The tests described in this chapter are in addition 
to the tests conducted to verify the accuracy of the programs. In all 
cases » phase corrections were performed. 

4.1 SOURCE INTENSITY MAP 

One of the important aspects of the plane wave tube is the behavior 
of the speaker array. It iS desirable for all speakers to produce 
identical outputs with the same phase angle. Also the spectrum produced 
by the speakers should be flat for a random white noise excitation. 

During the initial calibration tests of the test facility, it was 
concluded (Reference 2 ) that the incident wave can be considered plane 
only up to 800 Hz. With the acoustic intensity technique, this aspect 
can be easily verified. To determine the sound field characteristics 
of the test facility, an acoustic intensity survey was carried out 
along the cross section of the plane wave tube. The test facility has 
a cross section of 18 inches by 18 inches. Tests were conducted to 
measure intensity every two inches, using the procedures outlined in 
Appendix A. This gave intensity values of 81 grid points. During these 
tests, the gain values at the frequency ranges of the equalizer were 
set to zero. 
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The results of the tests are plotted in Figure 4.1 and 4.2, for 
300 Hz and 1000 Hz, respectively. The result?; are also available for 
every 1.25 Hz up to 500 Hz, and for every 12.3 Hz from 500 Hz up to 
5000 Hz. The software programs developed seem to work well for the 
type of analysis being done. From the tests, it was found that the 
number of grid points needs to be increased at high frequencies to 
obtain a good quality Intensity map. 

From Figure 4.1, it can be seen that two speakers (//2 and #6) 
are producing less power (10 dB lower than the other speakers). This 
phenomenon was seen at frequencies from 250 to 4000 Hz. Thereafter, 
these speakers behaved normally. But for these two areas, the result 
was reasonably flat. At 1000 Hz, the variations were much more severe. 
This could be due to the cavity resonances present in the test facility 
In general, the intensity was higher around the edges than at the 
center. The reason for this is not fully understood. However, based 
on this test, it is concluded that the KU-FRL acoustic test facility 
cannot be considered a plane wave facility above 1000 Hz. 

4.2 INTENSITY MAP WITH ALUMINUM PANEL 

At the KU-FEIt test facility, a 0.032" aluminum panel is used as 
the standard panel. The transmission loss (or noise reduction) values 
obtained with this panel are used for calibration. To determine the 
acoustic intensity characteristics of this panel, an intensity survey 
was carried out at the same 81 grid points as before, this time r;ith 
the 0,032" aluminum panel installed between the source and the 
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Source Intensity Map at 1000 Hz 



microphones. Figures 4.3 and 4.4 show the results at 300 Hz and 1000 
Hz, respectively. At 300 Hz, the intensity variation was within 10 dB 
at all points. At 1000 Hz, while the maximum variation was only 20 dB, 
the actual intensity value was 40 dB. It is anticipated that this low 
value of transmitted intensity may pose problems in accurate estimation 
of the intensity, especially if the panel exhibits higher transmission 
loss characteristics. This aspect was expected. At higher frequencies, 
the transmission loss will be higher because of the mass Jaw. Several 
methods are being studied to overcome this problem. They are installation 
of amplifier in the measurement channel, increasing the input signal 
strength, and finally filtering away the low frequency in the excitation 
signal using high-pass filters and then amplifying the signal. The 
third method will involve performing each test twice: once at low 

frequency, say up to 500 Hz; and the second time,, from* 500 Hz to 5000 Hz, 
This was how the tests were being done with the old measurement 
procedures. 

4.3 TRANSMISSION LOSS OF PANELS 

To compare the values of the measured transmission loss values 
obtained using this procedure with theoretical values, two pane 
were tested: a 0.032” aluminum panel and 40 oz/sq yd leaded vinyl. 

These specimens were tested at the KU-FRL acoustic test facility using 
the test procedures outlined in Appendix A. The resulting transmission 
loss characteristics are compared with the mass law. The behavior of 
the test panels is illustrated in Figures 4.5 and 4.6. 
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Figure 4.4; Intensity Map with 0.032" Aluminum Panel at 1000 Hz 




PANELi 0.032" ALUMINUM PANEL 
WEIGHTi 1.03 LB 
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Figure 4,5; Transmission Loss Characteristics of 0.032” Aluminum Panel 



original page 
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Figure 4,6: Transmission Loss Characteristics of 40 oz/sq yd Leaded Vinyl 



The transmission loss (TL) curve for a particular panel is obtained 
in the following manner. An intensity level survey is conducted without 
any panel. The intensity level is integrated over the entire panel area 
and the Incident sound power level is estimated. Then the tests are 
repeated with the panel installed between the source and the microphones 
and the transmitted sound power level is estimated. The difference in 
sound power level between these two measurements will give the 
transmission loss. This test procedure is similar to the measurement 
of insertion loss. However, because the intensity is measured close 
to the source and the .?anel is thin, the difference between this 
procedure and the two-room transmissicn loss measurement method is 
expected to be small. This procedure had to be adopted because the 
measurement of incident intensity alone on the source side is not 
possible with this technique. The reflected intensity from the test 
specimen will affect the Intensity being measured. This problem 
existed even in the old procedure. However, no correction was made 
for this; and the resulting measurements were termed noise reduction 
(and not transmission loss). 

Mass law values are also plotted in Figures 4.5 and 4.6, The 
mean square transmission loss measured in each case is approximately 
2 to 5 dB Iqwer at 5000 Hz. The integration of the several intensity 
values for both with and without panels, used in this procedure, is 
expected to yeidl an average transmission loss value. This is in 
comparison to the use of just one microphone situated close to both 
source and receiving side, which will result in position-dependent 
transmission loss values. Hence a difference between the two 


measurements was anticipated. The results with the old test procedure 
gave up to 6 dB higher noise reduction than the calculated noise reduction 
values. Hence, one-to-one comparison between the present procedure and 
the old single-microphone procedur^^ Is not considered valid. 

In addition to the lower-than-mass-law values observed at 5000 Hz, 
another aspect that will have to be studied Is the near-zero transmission 
loss observed around 300 Hz. The reasons for this are not fully 
understood and are currently being Investigated. The tests will be 
performed after modifications have been made to the MFD and the receiver 
cavity. These modifications include filling up the cavities with loose 
foam, moving wire gauze closer to the MPD, and Installation of the 
microphone in the MPD in such a way as to reduce the reflections 
from the support rod. Also, increasing the signal strength at high 
frequencies Is expected to Increase the dynamic range of measurements. 
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CHAPTER 5 


CONCLUSIONS AND RECOMMEKDATIONS 

In this report, the development , of a method to measure the 
transmission loss characteristics of panels using the acoustic Intensity 
technique is presented. This report also describes the theoretical 
formulation, installation of the hardware, modifications to the test 
facility, and development of computer programs and test procedures. 
Initial test results are also presented. Based on the tests conducted 
so far, the following conclusions have been reached. 

• The acoustic intensity technique can be adapted to measure 
the transmission loss characteristics of panels. Use of 
this method will give average transmission loss values as 
opposed to the position-dependent values obtained from 
single-microphone measurements. 

• The same technique and installation can also be used to plot 
the intensity maps of vibrating panels. Use of the microphone- 
positioning device greatly simplifies correct grid positioning. 

• The acoustic intensity programs can easily be written on a 
microcomputer, (Total cost of the microcomputer is less 
than $2500.) 

• The initial results indicate that transmission loss values 
measured using this method -are lower than theoretically 
predicted values. 

• This facility cannot be considered a plane wave facility at 
high frequencies above 800 Hz. 
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Based on the experience gained so far In using this technique, 
the following recommendations are made: 

• The computer programs should be optimized to lessen the test 
time and the memory storage requirements. This may require 
programming in assembly language or in higher level languages 
like "C." 

• The test facility should be modified to eliminate possible 
reflections from the microphone holding device. 

• The test procedure should be modified to increase the dynamic 
range at high frequency. This can be achieved by performing 
the tests twice, once at low frequency ('^500 Hz) and the 
second time at high frequency ('^'5000 Hz). During the high 
frequency testing, the low frequency component should be 

, eliminated by filtering. 

• At present, no spatial averaging is possible. The microphone 
positioning device can easily be modified to liave spatial 
averaging . 
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APPENDIX A 


ACOUSTIC INTENSITY TEST OPERATOR'S MANUAL 
A.l INTRODUCTION 

In this appendix the procedures for performing an acoustic intensity 
test are presented. All pertinent data are included for calibrating 
equipment, performing a test, and reducing the resulting data. However, 
at the time of this writing, the testing procedures are still being 
developed, and the following is not expected to be the final procedure. 
This manual will be updated as testing procedures are refined. 

The sections presented below are arranged in the order in which 
they have to be performed, although not all the sections will be required 
for every test. (For instance, the "infrequent” phase calibration tests 
are only required when daily calibrations using method /‘I do not agree 
with previous infrequent calibration test results.) The manual presents 
the procedures in a checklist format and does not dwell on why certain 
steps are required, as this would require an intimate knot^ledge of the 
system (which at this time is not fully developed). The manual is 
intended to serve as a quick reference guide during testing. 

A. 2 INFREQUENT MICROPHONE PHASE CALIBRATION PROCEDURE 

This section describes the procedures for performing an "infrequent" 
calibration test of the intensity test set-up for purposes of statistical 
assurance. That is, the daily calibration tests are compared to the 
results of an infrequent’ calibration test to determine the statistical 



validity of the daily calibrations. To do this, a large , random sample 
of calibration tests is required to estimate the calibration values of 
the population. The sample size required by the STAT.BAS statistics 
program is thirty tests. It is therefore necessary to generate thirty 
calibration tests as described below. 

1. Connect the 2->lnch calibration speaker to the metal calibration 
test tube with masking tape. 

2. Remove the red speaker wires (2) from the speaker panel and 
attach them to the black wires on the 2-inch speaker. 

3. Insert microphone //I into slot #1 in the wooden microphone 
support block, and microphone j/2 into slot #2 in the block. 

4. Insert microphone #1 into slot #1 in the Plexiglas microphone 
support ring, and microphone //2 into slot j/2 in the ring. 

5. Insert the Plexiglas ring into the metal test tube. Line up 
the markings on the ring and the tube, and secure together 
with masking tape. 

6. Attach a cable to the binary noise-out (A) post on the FFT 
analyzer. 

7. Connect the cable from the FFT to the speaker switch cable. 

8. Connect the cable from output 1 on the microphone power supply 
to channel A on the FFT analyzer. 

9. Connect the cable from output 2 on the microphone power supply 
to channel B on the FFT analyzer. 

10. Boot the Z-100 with the system/ testing disk in drive A, and 
a working disk in B. 

11. Turn on the following equipment: 
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a 


660B, 



b. mike power supply, 

c. equalizer, 

d. amplifier. 

Allow equipment to warm up for approximately 10 minutes. 

12. While the system is warming up, do the following: 

a. On Z-100: 

1) Type ZBASIC <ret>, 

2) Press F3 and when ’LOAD"’ appears, type CAL660 <ret>, 

3) When ’OK’ appears on the screen, type RUN <ret>. 

The computer will then ask for a file name. Don’t answer 
this question yet. 

b. On 660B: Make sure that the following lights are on. 

1) mode: Dual IK, 

2) function: TF, 


3) 

average : 

SUM, 

4) 

display: 

XPROP, <j)MAG, 

5) 

capture control: CONT, 

6) 

amplitude 

units: V/ABS 


7) scales: AMPL-LOG. 

13. Set the number of averages using the following sequence: 

a. 1-2-8, 

b. SET /STORE, 

c. AVGN. 

Check to see if 128 AVGN appears on the screen of the 660B. 

If it does, press CLEAR and repeat step 13. You are now ready 
to begin testing. 
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lA. Oa the Z-100, type in the file name In the following format. 
Do not hit return. 


Denotes Drive B 


File Name 

JL, 

■B:CALPXX.D(LO or HI) 


•Test number (01 to 30) 


D for Data 


/ 7 t 

ata / HI for High Frequency Callbrat'lons 

TAtf T r> ' (Up to 5000. Hz) 

LO for Low Frequency Calibrations ^ 

(Up to 500. Hz) 


15 . On the 660B : 

a. Set desired frequency range (500 “ LO, 5000 = HI). 

b. Turn on the speakers by putting the speaker switch In 

the "SO" position. (Note: speakers may be turned off 

whenever data are not being taken.) 

c. Adjust the analyzer channel amplitudes to maximum 
without overloading. 

d. Press START. 

e. Press RECALL and then AVON. 

f. When the test is done, the START light will go out. 

Turn off the speakers by ratuming the speaker switch 
to OFF. 

16. Press RETURN on the Z-100. This will initiate the transfer 

of data from the 660B to the Z-100. The transfer Is completed 

when the Z-100 asks for a new file name. 

17. On the 660B, press RCL D so that the function light goes out. 

18. Begin the next test by following the steps from step 14 onward. 

19. When finished, turn off: 

a. microphone power supply, 

b. amplifier. 


i- ! 


» 4 
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c. equalizer, 

d. 660B. 

e 

20. On the Z-100, type E <ret>. They type SYSTEM <ret>. 

21. Turn off the Z-100 and its monitor. 

22. Remove disks. 

23. Disconnect all cables and reconnect the speaker wires to the 
speaker panel. 

A. 3 PROCEDURE FOR RUNNING STAT.BAS: A PROGRAM TO REDUCE PHASE 

CALIBRATION DATA USED FOR STATISTICAL ASSURANCE 

This section describes the operation of the statistical analysis 
program STAT.BAS. This program calculates the sample mean and a confidence 
interval about the mean of thirty tests. The calculations are done such 
that 'the probability of the population mean falling within the confidence 
interval is ^3 percent. The program is designed to operate on data 
files of arbitrary length up to 403 values (that is, thirty files, each 
up to 403 values in length) . 

The following is a description of the preparation required to run 
STAT.BAS. 

Thirty calibration tests are required for both the low frequency 
(0-500 Hz) and high frequency (0-5000 Hz) tests. All thirty of these 
tests for either high or low ranges will fit on one disk. File names 
must be sequential, although it is not necessary to have them arranged 
in order. File name format must be as described in the previous 
section. 
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Program STAT.BAS requires a catalog file containing the names of 
all the calibration tests. For the low frequency tests, this catalog 
file Is ‘FILCAT.LFR; and It must reside on the same disk as the low 
frequency calibration tests. For the high frequency tests, the catalog 
file name Is FILCAT.HFR; and It must reside on the high frequency test 
disk. 

FILCAT.LFR and FILCAT.HFR are most easily generated by running 
simple utility programs called CATl.BAS and CAT2.BAS. CAT1.B>\S generates 
30 low frequency file names within the FILCAT.LFR file. These take the 
form 

B:CALP01.DLO 

B:CALP02.DLO 

B;CALP03.DL0 

• 

B:CALP30.DLO 

Similar file names are generated under FILCAT:HFR for the high frequency 
tests. 

STAT.BAS will open a FILCAT file and read in the calibration file 
names. STAT.BAS will open each calibration file in turn and read in 
the data in each file. When finished, the results will be written to 
a user-defined file in the form 

• 

• 

1.093 0.023 

1.110 0.023 

1.023 0.019 

0.913 0.006 
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The first column represents the mean phase difference between the 
two microphones (mike ifl with respect to mike //2) . The second column 


represents the confidence interval (the Interval is the mean +/- 
confidence interval value). Each row corresponds to a different 
frequency. For low frequency tests, the difference between txro 
consecutive rows is 1.25 Hz, For the high frequency tests, the 
difference is 12.5 Hz. Four hundred two (402) values are stored, 
starting at 0 Hz. 

The following is a step-by-step description of how to run STAT.BAS. 

1. Boot Z-100 with the system/ testing disk in drive A, and 
the low frequency calibration test disk in B. 

2. Type ZBASIC <ret>. 

3. When Z-IC / ^’esponds with 'OK,’ type L0AD"CAT1 <ret>. 

4. When Z-100 responds with 'OK,' type RUN <ret>. 

5. Z-100 will print 30 file names. When it is finished, type 
NEW <ret>, 

6. Replace the disk in drive B with the disk containing the 
high frequency tests. 

7. Type L0AD''CAT2" <ret>. 

8. When Z-100 responds with 'OK,' type RUN <ret>. 

9. Z-100 will print 30 file names. When it is finished, type 
NEW <ret>. The FILCAT files have now been created on their 
respective disks. 

10. Type L0AD"STAT" <ret>. 

11. When Z-100 responds with 'OK,' type RUN <ret>. 

12. Z-100 will ask 'ENTER NAME OF INPUT FILE CATALOG?' 


Type B:FILCAT.HFR <ret> 


13. Z-100 will ask ’ENTER OUTPUT FILE NAME?' 

Type STATHI.DAT <ret>. 

» 

14. Z-100 will ask 'ENTER ANALYSIS RANGE?' 

Type 5000 <ret>. 

15. Z-100 will now be busy for approjciniately 20 minutes. As 
each file is opened, its name will be printed on the monitor. 
When finished, Z-100 will beep three times. 

16. When Z-100 is finished, replace the disk in drive B with 
the disk containing the low frequency calibration tests. 

Type RUN <ret>. 

17. The answers to the three Z-100 questions asked in steps 
12, 13, and 14 are, respectively, 

B:FILCAT.LFT 

STATL0.DAT 

500. 

18. When Z-100 is finished, type SYSTEM <ret>. 

19. The high and low frequency statistical data now reside on 
the disk in drive A under the file names STATHI.DAT and 
STATL0.DAT, respectively. 

20. Copy files STATHI.DAT and STATLO.DAT to Data Disk #1. 

Do this using the following steps. 

a. Insert Data Disk #1 in drive B. 

b. Type STATHI.DAT B;*.* <ret>. 

c. Type STATLODAT B:*.* <ret> when Z-lOO answers with 'A:'. 

21. Turn off Z-100 and its monitor. 

22. Remove all disks. 
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A. 4 DAILY MICROPHONE PHASE CALIBRATION PROCEDURE; METHOD I 


This section describes the procedures for performing a daily 

m 

calibration test of the tvo B & K microphones. The procedure should 
be followed in the order given. Two tests are required; a low 
frequency test, and a high frequency test. Note that this phase 
calibration procedure is one of two procedures that may be used. 

1. Connect the 2-inch speaker to the test tube with masking tape. 

2. Remove the red speaker wires (2) from the speaker panel and 
attach them to the black wires on the 2-inch speaker, 

3. Insert microphone #1 into slot //I in the wooden microphone 
support block, and microphone //2 into slot #2 in the block. 

4. Insert microphone //I into slot //I in the Plexiglas microphone 
support ring, and microphone #2 into slot //2 in the ring. 

5. Insert the Plexiglas ring into the metal test tube. Line up 
the markings on the Plexiglas ring with those on the tube. 
Secure the ring and tube together with masking tape, 

6. Attach a cable to the binary noise-out (A) post on the FFT 
analyzer. 

?. Connect the cable from the FFT to the speaker switch cable, 

8. Connect the cable from output 1 on the microphone power 
supply to channel A on the FFT analyzer. 

9. Connect the cable from output 2 on the microphone power supply 
to channel B on the FFT analyzer, 

10. Boot Z-100 with the system/ testing disk in drive A, and 
Data Disk //l in drive B. 
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11. Turn on the following equipment: 

a. 660B, 

b. mike power supply, 

c. equalizer, 

d. amplifier. 

Allow equipment to warm up for approximately 10 minutes. 

12. While the system is warming up, do the following: 

a. On Z-100: 

1) Type ZBASIC <ret>, 

2) Press F3 and when 'LOAD”' appears, type CAL660 <ret>. 

3) When 'OK' appears on the screen, type RUK <ret>. 

The computer will then ask for a file name. Don't answer 
this question yet. 

b. On 660 b : . Make* sure that the following lights are on. 


1) 

mode; Dual IK, 


2) 

function: TF, 


3). 

average: SOM, 


4) 

display: XPROP, 

<^MAG, 

5) 

capture controls 

CONT, 

6) 

amplitude units: 

V/ABS, 


7) scales: AMPL-LOG. 

13. Set the number of averages using the following sequence: 

I * 

a. 6-4-0, 

♦ 

b. SET/STORE, ^ 

c. AVGN. 

'T 

Check to see If 640 AVGN- appears on the scrreen of the 660B. 
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If It does not, press CLEAR and repeat step 13. You are 
now ready to begin testing. 

14. On the Z-100: Type in the file name in the following format. 

Do not hit return. 


Denotes Drive B 


File Name (Transfer Function Calibration) 




:TFNCAL.D(LO or HI) 


f 

D for Data | HI for High Frequency Calibrations 

T_ . T r- (Up to 5000. Hz) 

LO for Low Frequency 

Calibrations (Up to 500. Hz) 


I 


15. On the 660B: 

a. Set desired frequency range (500 =» LO, 5000 = HI). 

b. Turn on the speakers by putting the speaker switch in 

the "SO" position. (Note: Speakers may be turned off 

whenever data are not being taken.) 

c. Adjust the analyzer channel amplitudes to maximum without 
overloading. 

d. il'ress START. 

e. Press RECALL and then AVON. 

f. When the display reads 128 AVON, press STOP. 

g. At this point, disassemble the microphone/block 
combination and repeat steps 3 to 5. 

h. After completing 15a to 15 g, press CONT on the 660B, 

i. When the display reads 256 AVON, press STOP, and repeat 
15g to 15h. Do this for 384 AVGN and 512 AVON also. 

j. When the test is done, the START light will go out. 

Turn off the speakers by returning the speaker switch 
to OFF. 


61 


16. Press RETURN on the Z-100. The transfer is completed when 
the Z-100 asks for a new file name, 

17. On the 660B, press RCL D so that the function light goes out. 

18. Begin the next test by following the steps from step 14 onward. 

19. When finished, turn off: 

a. microphone power supply, 

b. amplifier, 

c. equalizer, 

d. 660B. 

20. On the Z-100, type E <ret>. Then type SYSTEM <ret>. 

21. Type CALII <ret>. When Z-100 asks for the frequency 
range, type 500. <ret>. 

22. Z-100 will now do a statistical comparison of the phase 
values determined in the daily calibration with those from 
the Infrequent calibration tests. Whenever a value is found 
to be outside accepted limits, the Z-100 will give the 
frequency value and the value of the test statistic. Accept 
the calibration value with a <ret> if one of the following 
conditions is satisfied: 

a. Frequency value is less than 60 Hz or greater than 4000 Hz, 

b. The test statistic is close to 2.0 (within, say, five 
percent), 

c. iiie test statistic is less than 2.5 for less than five 
consecutuve frequency values, 

d. The test statistic is considerably greater than 2.0 but 
is an isolated excursion; that is, a single event, rather 
than a string of excursions greater than 2.0. 
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If the test statistic fails, then the calibration will have 
to be repeated. If a satisfactory calibration test cannot 
be made within three or four trials, then Section 
should be redone. 

23. If the values are acceptable, repeat steps ^^1 and 22 but for 
frequency value 5000. Hz. 

A. 5 MICROPHONE CALIBRATION USING SWITCHING TECHNIQUE TEST PROCEDURE 

CALIBRATION METHOD II 

This section describes the procedures for microphone channel 
calibration using the switching technique. 

1., Place the microphone positioning device (MPD) between the 
Beranek tube and the speaker wall. 

2. Insert the two microphones through the hole in the top of 
the MPD . 

3. Clamp the microphones into their respective holders (//2 mike 
closest to the speakers) mth a one-inch spacing for high 
frequency tests, and a four-inch spacing for low frequency 
tests. The microphones should be inserted so that they are 
upright ill their holder. 

4. Position the mikes in the upper left hand side (when facing 
the speakers) of the MPD, somewhere in the center of the 
cross section. The exact location is not important; but 
once it is in place, the MPD should not be moved. 

5. Place putty into the hole in the MPD around the microphone 


cables 



6. Leave the panel space empty. Close the tube and secure the 
MPD In place with the four "hook" bolts. 

7o Attach a cable to the binary noise-out (A) post on the FFT 
analyzer. 

8. Connect the cable from the FFT to the speaker switch cable. 

9. Connect the cable from output 1 on the microphone power supply 
to channel A on the FFT analyzer. 

10. Connect the cable from output 2 on the microphone power supply 
to channel B on the FFT analyzer. 

11. Boot Z-100 with the system/testing disk in drive A, and Data 
Disk #1 in drive B. 

12. Turn on the following equipment! 

a. 660B, 

b. mike power supply, 

c. equalizer, 

d. amplifier. 

Allow equipment to warm up for approximately 10 minutes. 

13. While the system is warming up, do the following: 

a. On Z-100: 

1) Type ZBASIC <ret>, 

2) Hit F3 and when 'LOAD"* appears, type INT660 <ret>, 

3) When 'OK' appears on the screen, type RUN <ret>, 

4) The computer vri.ll then ask for a file name. Do not 

answer this question yet. 

b. On 660B: Make sure that the following lights are on. 
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1) 

mode; DUAL IK, 


2) 

function; XSPEC, 


3) 

average; SUM, 


4) 

display; XPROP, 

(j)MAG, 

5) 

capture control; 

CONT, 

6) 

amplitude units; 

V/ABS, 

7) 

scales; AMPL-LOG 

• 


14. Punch the following buttons on the 660B in the order given; 

a. 2-5-6, 

b. SET/ STORE, 

c. AVON. 

Check and see if 256 AVON appears on the screen of the 660B. 
If it does not, hit CLEAR and redo the above sequence. 

You are now ready to begin testing. 

15. On the Z-lOO; Type in the file name in the following format. 
Do not hit return. 


Denotes Drive B 

B:XSNORM.D(LO or HI) 

t 

XS for Cross-Spectrum Low or High Frequency Calibration 

(LO = 500. Hz, HI = 5000. Hz) 

16. On the 660B: 

a. Set desired frequency range (500 or 5000 Hz). 

b. Turn on the speakers by putting the speaker svd.tch in 
the ”S0** position. 

c. Adjust the analyzer channel amplitudes to maximum without 
overloading. 
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d. Press START. 

e. Press RECALL and Chen AVGN. 

£. Note In Che Cesclng log Che InpuC chaiiuel maximum 
ampllCude seCClngs, microphone spacing, number o£ 
averages. 

17. When Che CesC is done, che START lighC will go ouc. Turn 
off Che speakers by reCurning Che speaker swiCch Co OFF. 

18. Press RETURN on Che Z-100. This will iniCiaCe Che Cransfer 
of daCa from 660B Co che Z>-100. The cransfer is compleced 
when Che Z-100 asks for a new file name. 

19. On che 660B, press RCL D so ChaC Che funcCion lighC goes ouC. 

20. Unscrew Che 4 holding bolCs. WiChouC changing Che posicion 
of Che MPD, swiCch Che locaCion of Che microphones. AfCer 
swiCching, f/1 microphone should be closer Co Che source. 

21. On Che Z-100, cype file name: 

B:X3SWCH.D(L0 or HI) 

Do noC hie reCurn. RepeaC seeps 16-19. 

22. RepeaC CesC for Che ocher frequency range (sCeps 16-21). 

23. Type E <rec> and SYSTEM <reC> on Che Z-100, and Curn off che 
following equipmenC : 

a. microphone power supply, 

b. amplifier, 

c. equalizer, 

d. 660B. 

24. Type INTCAL <rec>. When Z-100 asks for Che frequency range, 
Cype 500. <rec>. 
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25. Z-100 will now calculate the phase correction values and 
store the result in B:INTCAL.D(LO or HI). 

26. Repeat step 24 for frequency range value of 5000. Hz. 

A. 6 DAILY MICROPHONE MAGNITUDE CALIBRATION TEST PROCEDURE 

This section describes the procedures for perfonaing a daily 
magnitude calibration test of the two B & K microphones. The steps 
should be followed in the order given. Special note: Microphone //I 

corresponds to analyzer channel A, and microphone //2 corresponds to 
analyzer channel B. 

1. Remove the microphone connected to channel A from the MPD. 

2. Place the microphone into the connecting ring of the 
•'Pistonphone . 

3. Boot Z-100 with the system/testing disk in drive A, and 
Data Disk //I in drive B. 

4. Turn on the 660B. 

5. Turn on the microphone power supply. 

6. Allow the system to warm up for approximately 10 minutes, 

7. IThile the system is warming up, do the following: 
a. On the Z-100: 

1) Type ZBASIC <ret>. 

2) Press F3, and when 'LOAD'" appears, type PSP660 <ret> 

3) When *0K' appears, type RUN, 

The computer will then ask for a file name. Do not answer 
this question yet. 

*Briiel & Kjaer, Type 4220 
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b. On 660B: Make sure that the £ollovd.ng lights are on. 

1) mode: DUAL IK, 

b) function; PWR SPEC, 

c) average: SUM, 

d) display: MAG (either A or B) , 

e) capture control; CONT, 

f) amplltued units: V,ABS, 

g) scales; AMPL-LOG. 

8. Set the number of averages using the following sequence: 

a. 2-5-6, 

b. SET/ STORE, 

c. AVGN. 

Check and see if 256 AVGN appears on the screen of the 660B. 

If it does not, press CLEAR and redo the above sequence. 

9. You are now ready to begin testing. Do the following: 

a. On the Z-100: The Z-100 will ask ’FILE NAME FOR 

CHANNEL A?’. Type in the file name A: CALA.DAT. 

Do not enter a return yet. 

b. On 660B; Set the frequency range to 500 and press RECALL 
and then AVGN. 

10. Turn on the Pistonphone. 

11. Hit the START button on the 660B. 

12. When the test is done, the START light will go out. When 
this happens, turn off the Pistonphone. 

13. Hit RETURN on the Z-100. 

14. The transfer is complete when the Z-100 asks for a new file name. 
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15. On the 660B, hit RCL D so that the light goes out. 

16. Replace microphone A with microphone B in the Pistonphone, and 
output from microphone B should be connected to channel B 

on 660B. 

17. Begin the next test by following the steps from step 9a onward. 
This time the file name will be A:CALB.DAT. 

18. Type SYSTEM <ret>. When A; appears, type MAGCAL <ret>. 

19. The computer will then ask a series of questions. Answer 
each of them with the appropriate answer. 

20. When A: appears on the screen, type COPY A:CALDAT.DAT BiCALDAT.DAT 
<ret>. 

21. Turn off: 

a. microphone power supply, 

b. 660 b, 

c. Z-100 and its monitor. 

22. Remove the disks. 

23. Place the microphone back in the MPD and put away the 
Pistonphone. 

A. 7 ACOUSTIC INTENSITY TESTING PROCEDURE 

This section describes’ the procedures for running an intensity test. 
Follow these steps in the order given. Note: Each disk will hold only 

35 files. You will need several data disks for these tests. 

1. Place the microphone positioning device (MPD) between the 
Beranek tube and the speaker wall. 
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2. Insert :he two microphones through the hole In the top of the 
MPD. 

3. Clamp the microphones into their respective holders 02 mike 
closest to the speakers) with a one-inch spacing for high 
frequency tests, and a four-inch spacing for low frequency 
tests. The microphones should be Inserted so that they are 
upright in their holder. 

4. Position the mikes in the upper left-hand side (when facing 
the speakers) of the MPD, one inch from the left side and top 
side of the inner surfaces. This represents position (1,1). 
While holding the positioning cranks fixed, move the position 
indicators (springs) so that they are centered over the 
indicator scales at the value (1,1). 

5. ' Place putty into the hole in the MPD around. the microphone 

cables. 

6. If a "receiver side" test is to be done, place the panel to 
be tested in the MPD. If a "source side" test is to be done, 
leave the panel space empty. Close the tube and secure the 
MPD in place with the four "hook" bolts. 

7. Attach a cable to the binary noise-out (A) post on the FFT 
analyzer. 

8. Connect the cable from the FFT to the speaker switch cable. 

9. Connect the cable from output 1 on the microphone power 
supply to channel A on the FFT analyzer. 

10. Connect the cable from output 2 on the microphone power supply 
to channel B on the FFT analyzer. 
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11. Boot Z-100 with the system/ testing dish in drive A, and 
an empty working disk in B. 

12. Turn on the following equipment; 

a. 660B, 

b. mike power supply, 

c. equalizer, 

d. amplifier. 

Allow equipment to warm up for approximately 10 minutes. 

13. While the system is warming up, do the following: 

a. On the Z-100: 

1) Type ZBASIC <ret>; 

2) Hit F3 and when ’LOAD*'* appears, type INT660 <ret> 

3) When '0K“ appears on the screen, type RUN <ret>, 

4) The computer will then ask for a file name. Do 
not answer this question yet, 

b. On 660B; Make sure that the following lights are on. 

1) mode: DUAL IK, 

2) function; XSPEC, 

3) average: SUM, 

4) display; XPROP, <{>MAG, 

5) capture control: CONT, 

6) amplitude units; V/ABS, 

7) scales: AMPL-LOG, 

14. Punch the following buttons on the 660B in the order given: 

a. 2-5-6, 

.b. SET/STOae,^ 

c. AVGN. 
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Check and see if 236 AVGN appears on the screen of the 660B. 
If it does not, hit CLEAR and redo the above sequences. 

You are now ready to begin testing. 

15. On the Z-100: Type in the file name in the following format. 

Do not hit return. 


Denotes Drive B ^Source or Receiver Side 

i 

B:XS(S or R)XX(L or H).XXX.< Test number 

t t 

XS for Cross-Spectrum 1 low or High Frequency Tost 

I (L - 500 Hz, H - 5000 Hz) 

Mike Position 

Examples: B;XSS31L.012 or B:XSS78H.005 

16. Begin the test by moving the cranks on the MPD so that it 

reads (1,1). The tests will need to be done in the following 
order: (1,1) - (1»9), (2,1) - (2,9), etc. After test (6,9), 

consult the special instructions listed at the end of this 
section. 

17. On the 660B: 

a. Set desired frequency range. 

b. Turn on the speakers by putting the speaker switch 
in the ”SO*' position. 

c. Adjust the analyzer channel amplitudes to maximum 
without overloading. 

d. Press START. 

e. Press RECALL and then AVGN. 

f. For each' position, note in the testing log the input 
channel maximum amplitued settings, microphone spacing, 
number of averages, and the number of the disk used. 
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18. When the test is done, the START light will go out. Turn off 
the speakers by returning the speaker switch to OFF, 

19. Press RETURN on the Z-100. This will initiate the transfer 
of data from the 660B to the Z-100, The transfer is completed 
when the Z-100 asks for a new file name. 

20. On the 660B, press RCL D so that the function light goes out. 

21. Begin the next test by repositioning the microphones to the 
next testing location and repeating steps 15 to 21. 

Remember that upon reaching position (6,9), the special 
instructions must be followed. 

22. After reaching position (9,9), and if no more testing is to 
be done, type E <ret> and SYSTEM <ret> on the Z-100, and 
turn off the following inuip’^Jent ; 

a. microphone power supply, 

b. amplifier, 

c. equalizer, 

d. 660B, 

e. Z-100 and its monitor. (Remove all disks.) 

*** Special Instructions *** 

These instructions are necessary at the time the microphones are 
moved past position (6,9). Beyond this location the preamp cables in 
the MPD begin to interfere with the operation of the device, and the 
microphones themselvas obstruct further progress. To avoid this, the 
microphones are turned over so that they face down. This operation is 


described below; 



1. Turn off the microphone power supply. 

2. Remove the "hook" bolts from the MPD and slide the tube back, 

3. Remove the aluminum beam to which the microphones are fastened 
by removing the bolt that fixes it to the Lucite mount. 

4. Turn the mikes face down and reattach the aluminum beam to 
its mount with the bolt. Make sure the arrangement of the 
mike cables allows for freedom of movement of the mikes. 

5. Position the microphones in the (9,1) position by carefully 
measuring one inch from the bottom and left side.* 

6. While keeping the positioning cranks fixed, move the position 
indicators (springs) so that they are centered over (9,,1) on 
the scales. 

7. Crank the mikes into the (7,1) position. 

8. Close the tube and bolt it together with the four "hook" bolts. 

9. Turn on the mike power supply. 

10. Proceed with testing from step 15. 
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LISTING OF COMPUTER PROGRAMS 


B.l LISTING OF PSP660.BAS 

ORIGINAL PAGE W 

5 REH PSPAAO OF POOR QUALITY 

10 SCREEN 0»0 
20 BEFINT I-N 
30 KEY OFRCLStCLOSE 
40 LOCATE 25 »1 

50 CI**ZZZ0A31SY251SY261SY272SyAJ6FlE0*8* 

60 PRINT STRING*(60f' •) 

70 SYN$=CHR4<22) 

80 LOCATE Ifl 
90 SPl',ED$**9600* 

100 CC«FIL$=*C0Mi:*+SPEEDi+'fEf7»l* 

110 OPEN COHFILi AS 41 

120 OPEN *SCRNf FOR OUTPUT AS #2 

130 LOCATE 25fi:PRINT '660B MAGNITUDE CAL TRANSFER PROGRAMS 

140 LOCATE IfltPRINT STRIN6»<60f* •) {LOCATE 1»1 

141 JC=1 

150 IF JC=1 THEN LINE INPUT 'FILE NAME FOR CHANNEL A? < E TO EXIT > {SDSKFILf 

151 IF JC=2 THEN LINE INPUT ‘FILE NAME FOR CHANNEL B? < E TO EXIT > tSDSKFILI 

152 IF JC»2 THEN MID$<C$»33)=*9' 

153 IF JC=2 THEN MIDi(C$f 29)=*2' 

160 IF DSKFIL$=‘E* THEN 500 

170 LOCATE IfllPRINT STRING$(60»* MILOCATE Irl 
180 OPEN DSKFIL$ FOR OUTPUT AS *3 

190 FIEL041f 3 AS Hl$»6 AS Fl$f3 AS El»»3 AS H2»f6 AS F2$*3 AS E2i>3 AS H3it6 
AS F3$f3 AS E3$f3 AS H4$r6 AS F4tf3 AS E4$f3 AS HS«f6 AS F5«f3 AS E5«i3 AS 
H6$r6 AS F6ir3 AS E64 
200 FOR IC=1 TO 33 
210 D$=HID4(C4rIC>l) 

220. GOSUB 340 
222 FOR IK=1 TO 75 
224 NEXT IK 
230 NEXT IC 
240 REM CONTINUE 
250 LOCATE Ifl 
260 GOSUB 370 
320 CL0SE*3IJC=JC+1 

330 IF JC =2 THEN GOTO 150 ELSE GOTO 500 

340 PRINT ♦1»D$J 

350 IF LOC<1)=0 THEN 350 

360 At=’INPUT$(lf4i){IF ASC(A$)«6 THEN RETURN ELSE PRINT #2» 'ERROR SENDING DATA* 

{STOP 

370 J«0 

380 FOR IC-1 TO 86 
390 PRINT ♦lfSYN$f 
400 J=J+1 

410 IF L0C(1)<72 THEN 410 
420 GET 41f72 

430 A*=Fl$+‘EnEl$+'f'+F2*+'E'+E24+'»'+F3l+'E'+E3*+'»'+F4t+'E*+E4$+'»'+F5$+'E' 
+E5»+'f'+F6$+'E"+E6$ 

440 IF J>67 THEN GOTO 460 
450 PRINT 43fA* 

460 NEXT IC 

470 IF L0C(1)<>2 THEN 470 
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B.2 LISTING OF MAGCAL.FOR 

ORIGINAL PAGE IS 
OF POOR QUALITY 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C HA6N1TUDE CALIBRATION PR06RAH 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c 

c 

Cf *«*«**«##♦« PR06RAHHER I R »NAVANEETHAN 
C««t****n*4* DATE t SEPTfl 83 

c«4*****Mm VERSION t 1 
c 
c 

C234567 

IMPLICIT REAL(Q-ZJ 
IMPLICIT REAL(A-F) 

IMPLICIT INTEGERS <I-N) 

DIMENSION CHA(402)fCHB(402) 

CHARACTER t\ HAfHB 
DATA HA/'AV 
DATA HB/*BV 

IF (I0REA0<6»2f0»*AICALA.DAT*)> 60 TO 10 
IF (I0READ(7»2f0f ‘AtCALB.DAT*)) GO TO 11 
J»l 

DO 1 1=1 r67 

READ (6fl00) (CHA<J1)» Jl=JfJfS) 

READ <7»100) (CHB(Jl)f Jl=JiJ45> 

J-J+A 

1 CONTINUE 

100 F0RMAT<6(E0»0)) 

IF (I0CL0S<6)> STOP 
IF (I0CL0S(7)) STOP 
YMCHA = CHA(l) 

YMCHB = CHBd) 

DO 2 I=5i400 

IF (YMCHA*6E.CHA(I)> GOTO 3 
YMCHA=CHA<!) 

IAMAX=I 

3 CONTINUE 

IF (YHCHB.6E.CHB<I>) GOTO 2 
YMCHB=CHB(I) 

IBMAX=I 

2 CONTINUE 
YA=0. 

DO 4 I=IAHAX-3fIAMAXP3 

4 YA=YA+CHA(I) 

YB=0. 

DO 5 I=IBMAX-3»IBMAX+3 

5 YB=YB+CHB(I) 

PCAL=1 0 ,*»< 1 24 . /20 00002 
AKCHA=PCAL/SQRT(YA) 

AKCHB=PCaL/SQRT<YB) 
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o o r> 


ORIGINAL PAGE IS 

URITEai200) OF POOR QUALITY 

200 FORMAT C' '» 'CHANNEL A DETAILS') 

HRITE(l»201) 

201 FORMAT (' '» 'ENTER MAX AMPLITUDE SETTING i '%) 

READ(lf300> AMPCHA 

300 FORMAT (F0*0> 

WRITE(lf202) 

202 F0RMAT(' 'i 'ENTER * OF AVERAGES t '♦> 

READ(lf301) NAV6A 

301 FORMAT(IO) 

«RITE(lf203) 

203 FORMAT(' CHANNEL B DETAILS') 

HRITE(lf201) 

READ(li300) AMPCH6 
WRITE(lf202> 

READ<lf301) NAVGB 

IF (lOURIT (8»2»0»*CALDAT.DAT’)> STOP 
URITE (8»205) AKCHAiAKCHB 

205 FORMAT(' '.2E15.5) 

WRITE(a»205) AMPCHAfAHPCHB 
URITE<8>206) NAVGAfNAVGB 

206 FORMAT<' 'r2I5) 

IF <I0CL0S<8>) STOP 
GOTO 6 

10 «RITE(1»500) HA 

500 FORMAT(' 'ERROR OPENING DATA FILE OF CHANNEL '>AD 
60 TO 6 

11 URITE(lf500) HB 

6 CONTINUE 

STOP 

END 
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B.3 LISTING OF TFN660.BAS 


t 

ORlGli'^AL PAGE 59 

5 REM TFN660.BAS OP POOR QUALITY 

10 SCREEN OfO 
20 DEFINT I-N 
30 KEY QFFJCLStCLOSE 
40 LOCATE 25^1 

50 C4»*ZZZ0A31SY251SY261SY272SUAJ3EOE4ZZZ*8* 

60 PRINT STRING$(60»' •> 

70 SYNi«CHRi<22) 

80 LOCATE 111 
90 SPEED4=*?600* 

100 C0MFIL*»*C0Hi:*+SPEED$+SE»7»l* 

no OPEN C0HFIL4 AS 41 

120 OPEN 'SCRMt* FOR OUTPUT AS *2 

130 LOCATE 25»nPRlNT *660B TRANSFER PROGRAMS 

140 LOCATE 1»1{PRINT STRINGi(60,* StLOCATE Ifl 

150 LINE INPUT 'INPUT FILE? <TYPE E TO EXIT > ‘.SDSKFIL* 

160 IF DSKFlLis'E* THEN 500 

170 LOCATE IfliPRINT STRING$(60r* SJLOCATE Ifl 
180 OPEN DSKFILi FOR OUTPUT AS 43 

190 FIELD4if 3 AS Hl«f6 AS F14f3 AS El«f3 AS H2$f6 AS F24t3 AS E2«f3 AS H34f6 
AS F34f3 AS E3$t3 AS H44f6 AS F44f3 AS E4$f3 AS H5«f6 AS F5«f3 AS E54r3 AS 
H6if6 AS F6if3 AS E6$ 

200 FOR ICsl TO 36 
210 D4-MID4(C4»ICfl) 

220 GOSUB 340 

222 FOR IK=1 TO 75 

224. NEXT 'IK 

230 NEXT IC 

240 REM CONTINUE 

250 LOCATE Ifl 

'Uts eesuB 370 

}l165*Z»{GbSUB 340 
%mi GOSUB 340 

tn !*=• = * JGOSUB 340 
TOO B$=*9'{G0SUB 340 
310 GOSUB 370 
320 CLOSE 43JCLS 
330 GOTO 130 
340 PRINT 4lfD4i 
350 IF LOC<1)=0 THEN 350 

360 A$<INPUTi(lf41)iIF ASC(A4)>6 THEN RETURN ELSE PRINT 42f 'ERROR SENDING DATA* 

:STOP 

370 J-0 

380 FOR IC=1 TO 86 
390 PRINT 41fSYN$f 
400 J=J+1 

410 IF L0C<1)<72 THEN 410 
420 GET 41f72 

430 A4*Fl$+*E*+El$+*f*+F24+*E*+E2$+*f*+F3$+*E*+E34+*f*+F44+*E*+E4$+*f*+F5$+*E* 
+E5*+*f *+F6$fE*+E6$ 

440 IF J>67 THEN GOTO 460 
450 PRINT 43fA$ 

460 NEXT IC 
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B.4 LISTING OF STAT.BAS 


original page m 

OF POOR QUALITY 


10 CL3 

20 OPTION BASE 1 
30 DIN X(513)fX2<513)rAVHA6<S13}r'v'AR<S13>fTAM(31)»CI(513) 

40 FOR 1=1 TO 31 
50 READ TAVKI) 

60 NEXT I 

70 DATA 12.71»4<303f3»182f2.774»2.S71>2.447»2.36St2.306 
00 DATA 2.262»2.228p2.201r2*179f2.160f2.145f2.131f2.12f2»llr2.101 
90 DATA 2*093>2.086t2.08f2.074»2.069r2»064>2.06»2.056f2*0S2>2.048 
100 DATA 2.04S»2»042»2.02 

no N=0 

120 INPUT 'ENTER NAHE OF INPUT FILE CATALOG '»A$ 

130 OPEN 'I'»*1»A$ 

140 INPUT 'ENTER OUTPUT FILE NAHE 'IDS 

150 INPUT 'ENTER SPECTRAL LINE SPACING SLS 

160 IF EOF(l) THEN GOTO 390 

170 N=N+1 

180 INPUTtlf N4 

190 PRINT N4 

200 OPEN 'I'f42» N$ 

210 INPUT#lf H$ 

220 PRINT M4 

230 OPEN 'I'i*3» H4 

240 1=0 

250 IF E0F(2) THEN GOTO 360 
260 1=1+1 
270 INPUT*2fA 
280 INPUT*3»B 

290 DEG=ATN(B/A)*<180/3»14159) 

300 IF(A<0)AND<B>0) THEN DEG=DEG+130 
310 IF(A<0)AND(B<0) THEN DEG=DEG-180 
320 X(I)=X(I)+DEG 
330 X2(I)=X2(I)+DEG"2 
340 PRINT ItLS-LSJDEGfX(I)fX2(I) 

350 GOTO 250 
360 CLOSE 43 
370 CLOSE 42 
380 GOTO 160 
390 CLOSE 41 
400 OPEN '0'f41»B$ 

410 K=N 

420 IF M>31 I'HEN K=31 
430 FOR J=1 TO I 

440 VAR<J) = (N»X2(J)-X<J)''2)/(N*(N-1)) 

450 X(J)=X(J)/N 

460 CI(J>=TAV(K)»(VAR(J)/N)"2 

470 PRINT J*LS-LSJX(J)JVAR<J)»CI(J) 

480 PRlNT41rUSIN6*44444*444 ' $ J»LS-LS*X( J) fCK J) 

490 NEXT 
500 CLOSE 41 
510 FOR 1=1 TO 3 


540 MEXT 
550 NEXT 
540 ENZi 


ORfGINAl PAGE 
OF POOR QUALITY 


B.5 LISTING OF CALII.FOR 


ORIGINAL PAGS I® 

OF POOR QUALITY 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C PHASE CALIBRATION PR06RAH FOR INTENSITY 

C TRANSFER FUNCTION NETHOB 

C TYPE I V 

C 
C 

C VERSION : 1 

C PR06RAHHER { R.NAVANEETHAN 

C DATE : 20-HAY-83 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C23AS67 

c 

CM**#*m CHANGE DIMENSION STATEMENTS DEPENDING UPON THE ANALYSER 
C 

COMPLEX CALA (S12) 

COMPLEX CMPLX 

CHARACTER «15 NAHETf NAMEl »HAHE2»NAME3rNAME4rNAMEC 
CHARACTER U AR 
REAL X<512)»Y(512) 

REAL ABIB(4) 

DATA ABIB/4«0./ 

DATA NAMEl/'BtTFNCAL.DLOV 
DATA NAHE2/*BtTFNCAL.DHIV 
DATA NAHE3/*B1CAL1I.DL0‘/ 

DATA NAME4/*BiCALII,DHIV 
RAD=180./3, 1415962 
URITEvlf600) 

600 FORMAT (' £NTER FREQUENCY RANGE t 
READ(lf601) SFREQ 

601 FORMAT(FO.O) 

IFLAG=1 

NAHET==NAME1 

NAMEC=NAME3 

IF < SFREQ *GT» 1000.) IFLA6»2 
IF (IFLAG.NE.2) GOTO 112 
NAMET-NAME2 
NAHEC=NAME4 
112 CONTINUE 
C 

C«0**44444* CHANGE N DEPENDING UPON THE ANALYSER 
C 

Ns402 

C 

cm*4**«*« OPEN DATA FILE CONTAINING TRANSFER FUNCTION DATA 
C 

IFa0REA0(8f2f0rNAHET>) STOP 
J*1 

DO 1 I=l»67 

READ(SflOO) (X(J1)» Jl»JfJ4S> 

JsJ+6 

1 CONTINUE 
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oooooooooo 


ORIGINAL PAG2 
OF POOR QUALITY 


J*1 

DO 11 1=1 >67 

READ(8»100) (Y<J1)> Jl»JfJ+5> 

J=Jf6 

11 CONTINUE 
C 

C*««******* CHANGE FORMAT STATEMENT DEPENDING UPON THE ANALYSER 
C 

100 FORMAT(6EO.O) 

IF(I0CL0S(8)) STOP 

namet='a;statlo»dat' 

IF (IFLAG*EQ*2> NAMET=*A:STATHI .OAT* 
IFa0READ(6f2>0>NAMET)) STOP 
READ(6fil4) FREQ 
114 FORHAT(FO.O) 

IF<FREQ.NE,3FREQ) GOTO 999 
BU=SFREQ»2. 56/1024. 

DO 3 I=1>N 

READ<6>110) AHEANfALUL 
110 F0RMAT(F9.3»lXfF9.3) 

IF(ALyL.LT.,l) ALVL=.l 
THETA1=ATAN2<Y(I)>X(I))»RAD 
IF((IFLA6.EQ.2) .AND.(I*6T.280)) GOTO 4 
IF(I.EQ.l) THETA1=0.0 
SIGMA=ALVUSQRT(30.) 

ZSTAT= ( AMEAN-THETAl ) / ( SIGMAtSQRT < 1 , /30 , + 1 , /5 . ) ) 

937 FORMATC' ' >F15.2f2F15.3) 

IF<ABS(ZSTAT).LE.1.96) GO TO 4 
«RITE(lr990) 

990 FORHAT(' FREQUENCY CAL VALUE ZSTAT 

FRQ=FL0ATa-l)*8U 
URITEa»987) FRQ»THETAlf ZSTAT 
989 URITE(1»113) 

113 FORMAT (' VALUES MOT WITHIN LIMITS! ACCEPT OR REJECT <A/R> 

READ<lr988) AR 
988 FORMAT (AO) 

IF<(AR.NE.*A*) ,AND.(AR.NE.*R*)) GOTO 989 
IF(AR.EQ.*R*) GO TO 998 
4 CONTINUE 

THE7A s*-THETAl/RAD 

CALAa)=CMPLX<COS(THETAl)»SIN(THETAl)) 

3 CONTINUE 

C 

C*******44# OPEN PHASE ''AL DATA FILE 
C 

IF(IOURIT(10»2»0»NAMEC)) STOP 

ABIB(1)=SFREQ 

ABIB<2)=2. 

WRITE(10»102)(ABIB(J)»J=1»4> 

102 FORMAT (' ^f4F15.5) 

DO 2 I=1»N 

URITE(10>101) CALA(I) 

101 FORMAT(' S2E15.5) 


2 CONTINUE 

IF<IOCLOSUO)) 5T0P 
GOTO 9999 

999 WRITEdr’OO) 

700 FORHAU' FREQUENCY MI8-HATCHM 
60 TO 9999 
998 CONTINUE 
9999 CONTINUE 
STOP 
END 


original page m 

OF POOR QUALITY 


I ■ 

V • 
^ . 





86 



B.6 LISTING OF INT660.BAS 


S REH INT660»BAS 
10 SCREEN OfO 


PAC3 Tif 
OF POOR QOAUW 


20 DEFIHT I-N 
30 KEY OFF ICLSt CLOSE 
40 LOCATE 25 »1 

SO C*=>'Z2Z0A31SY251Sy261SY272S«AJ2E0E4ZZZ=8* 

60 PRINT STRING$(60»* *) 

70 SYNI=CHR$(22) 

00 LOCATE 111 
90 SPEEn$=*9600' 

100 C0HFIL4=*C0M1{ •FSPEEDff iEi7f1 ' 

110 OPEN CONFILf AS ♦! 

120 OPEN *SCRN‘.‘ FOR OUTPUT AS #2 

130 LOCATE 25»ltPRINT *660B TRANSFER PROGRAHS 

140 locate 111 ‘.PRINT STRING*(60i* •) {LOCATE lil 

ISO LINE INPUT 'INPUT FILE? <TYPE E TO EXIT > {'iDSKFlLI 

160 IF DSKFILI='E* THEN 500 

170 LOCATE liltPRINT STRIN6$(60f* ') {LOCATE lil 
180 OPEN DSKFILI FOR OUTPUT AS *3 

190 FIELD#li 3 AS Hl$i6 AS F14i3 AS El»i3 AS H24i6 AS F2$i3 AS E2$i3 AS H3$i6 
AS F3$i3 AS E34i3 AS H4$i6 AS F4$i3 AS £4$f3 AS HS«i6 AS F5$i3 AS ES«i3 AS 
H6$i6 AS F6li3 AS E64 
200 FOR IC=1 TO 36 
210 Hi=HID$(C»iICil) 


220 QOSUB 340 
222 FOR IK=i TP 75 


224 NEXT IK 

230 NEXT IC 

240 REH CONTINUE 

2S0 LOCATE lil 

260 GOSUB 370 

270 DI='Z“,60SUB 340 

230 Di=*Z* {GOSUB 340 

290 D4='=' {GOSUB 340 

300 Dl='9* {GOSUB 340 

310 GOSUB 370 

320 CLOSE *3{CLS 

330 GOTO 130 

340 PRINT 41iD$i 

350 IF L0C(1)=0 THEN 350 

360 A$=INPUT$(li*l){IF ASC(A«)s6 THEN RETURN ELSE PRINT *2i'ERR0R SENDING DATA* 

{STOP 

370 J~0 

380 FOR IC=1 TO 86 
390 PRINT tliSYNSf 
400 J=J+1 

410 IF L0C(1K72 THEN 410 
420 GET ♦ii72 

430 A$=Fl$+*E'+El$+'i'+F24+'E'+E2$+'i'+F3l+'E*+E3$+'i'+F4$+'E'+E4$+'i'+F5i+'E' 
+E5$+'i'+F’6$+'E'+E6$ 


440 IF J>67 THEN GOTO 460 
450 PRINT 43iA$ 

460 NEXT IC 
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470 IF LOG (1)02 THEN 470 
480 A$=INPUn<2»*l> 

490 RETURN 

500 CLOSE {KEY ONtKEY OFF 


ORIGII\!AL PAGE JS 
OF POOR QUAUl'Y 
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B.7 LISTING OF INTCAL.FOR 


ORIGINAL 

OF POOR QUALITY 




I 

! 



ii ' 


ffW 

M 

U 


ii ^ 

U {t 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCC'CCCCCCCCCCCCCCCCCGCCCCCCCCCCCCCCCCCCCCCCC 


c c 

C PHASE CALIBRATION PROGRAM FOR INTENSITY C 

c c 

c c 

C VERSION : 1 c 

C PROGRAMMER i R.NAVANEETHAN C 

C DATE J 02-SEP-83 C 

C C 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C234567 


c 

ztn****** 

CC«***44t* 

Zt**H**U 

c 

c 

ctm**#** 

c*««**#*tt 

c 

c 

c 

c 

c#*#+***** 

c 

c 

cmm*** 

c 

c 

c 


THIS PROGRAM CALCULATES THE PHASE CALIBRATION USING 
THE METHOD DESCRIBED IN 'THE APPLICATION OF ACOUSTIC 
INTENSITY FOR ENGINE NOISE REDUCTION* BY H.D. CROKER ET AL 
PRESENTED AT THE INTERNATIONAL CONF. ON RECENT ADVANCES IN ACOU 
ACOUSTIC INTENSITY METHODS» SENLISi FRANCE^ OCT 1981 


IN THIS AND OTHER PROGRAMS ON ACOUSTIC INTENSITY PROGRAMS 
THE CROSS-SPECTRAL AND INTENSITY DATA ARE STORED IN TUO 
ARRAYS : 

A» DATA ARRAY IN UHICH DATA VALUES ARE STORED. 

B. IDENTIFICATION ARRAY IN UHICH THE SAMPLING 

FREDUENCYf DISTANCE BETUEEN MICS ETC ARE STORED. 


DIMENSION STATEMENTS 


DEPENDING UPON THE FFT ANALYSER CHARACTERISTICS THE 
DIMENSION VALUES NEED TO BE CHANGED. 


COMPLEX CALA(402)»CAL.B(402).C1 


C 

C#**** THE FOLLOWING STATEMENT IS PECULIAR TO SUPERSOFT FORTRAN AND 
Z*H** SHOULD BE OMITTED FOR OTHER FORTRAN COMPILER 
C 


COMPLEX CMPLX 
REAL X(402>»Y(402) 

REAL ABIB(4)»BBIB(4) 

CHARACTER *12 NAMElfNAME2» NAME3 
CHARACTER*2 STRl 

DATA namei/*b:xsnorm.dlo*/ 

DATA ;^!AMt2/*B:XSSUCH.DL0*/ 

DATA NAM£3/*B:INTCAL.DLQ*/ 

DATA STR1/*HI*/ 


c - ^ 
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ORIGINAL 

OF POOR QUALITY 

DATA ABIB/^CO./ 

DATA BBIB/4»0*/ 

C 

READ DATA FRQH X-SPECTRUH DATA FROM THE FFT ANALYZER 
C«4*#«m« STORED IN THE DISC 

Ct**U**** THIS FORMAT FOR OPENING THE DISC FILE IS PECULIAR 
Ct«*«««444 TO SUPER SOFT COMPILER 
C 

MRITE(lf700) 

700 FORMATr ENTER FREQUENCY RANGE OF ANALYSIS t '♦) 

READ (If 800) SFREQ 
800 FORMAT (F0»t>) 

IFLAG=1 

IF(SFREQ,LT»1001») GO TO 87 
IFLA6=2 

CALL PUTCHR < NAMEl » 1 1 f KH AR < S7R1 f 1 ) ) 

CALL PUTCHR ( NAME 1 f 1 2 f KHAR ( STR i f 2 ^ 

CALL PUTCHR ( NAME2 » 1 1 ? KHAR < STRl f 1 > ) 

CALL PUTCHR<NAHE2»12fKHAR(8TRlf2)) 

CALL PUTCHR < NAMES f 1 1 » KHAR ( STRl » 1 ) ) 

CALL PUTCHR (NAMES f 12 f KHAR< STRl » 2 ) ) 

37 CONTINUE 

IF aOREAD(6f2iOfNAHEi)) STOP 
C 

C«*44««*44 CHANGE N DEPENDING UPON THE ANALYSER 

C 

N=402 

J=1 

DO 1 1=1 r67 

READUflOO) (X(Jl)f J1=J»JFS) • 

JsJ+6 

1 CONTINUE 

J=»l 

DO 10 I=lf67 

READ(6flOO) (Y(Jl)f Jl-JfJ+5) 

J=J+6 

10 CONTINUE 

IF a0CL0S(6)) STOP 
C 

C«444t«A4# THE FOLLOWING FORMAT IS PECULIAR TO SUPER SOFT 
C****««44# FOR FREE FORMAT INPUTS SEPARATED BY COMMAS. 

C 

ICO FORHAT<AEO.O> 

C 

DO 12 I=lf402 
Y<I)« -Y(I) 

CALA(I) = CHPLX(X(I)fY(I)) 

12 CONTINUE 
C 

IF<I0READ<9i2f0fNAHE2)) STOP 
J*1 

DO 2 I=lf67 

READ (9fl00) (X<Jl)f Jl=J>JfS) 



r. I 

t t 
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J«J+6 POo»* 

2 CONTINUE 
J*l 

DO 20 I=l»67 

READ (9»100) <Y<J1)» J1=J»J^5) 

J=J+6 

20 CONTINUE 
IF(I0CL0S<9)> STOP 
DO 3 I=lfN 

3 CALB<I)=CHPLX(X(I)»Y(I)> 

C 

C«4**«* END OF DATA READ 
C 

IF(I0URIT(10r2»0fNAHE3)> STOP 
ABIB(i):=SFREQ 
BBIB(i)=SFREQ 
ABIB(2)=2. 

BBIB(2)=2, 

URITE(10fl02)(ABIB(J)>J=l,4) 

102 FORHAT(' 'f4El5»5) 

C 

CM**4« FOR MORE DETAILS SEE REF ABOVE 
C 

DO 21 I=1>N 
C 

C****mt# SEE REPORT KU-FRL-417-22 FOR DETAILS OF THIS REPORT. 
C 

Cl a CALA(I)/CALB(I) 

THETA = ATAN2(AIMAG(Cl)/REAL<Cl))/2. 
CALA(I)=CMPLX(COS(THETA)»SIN(THETA)) 

C 

C»*#**#«* WRITE TO DISC NEW PHASE CAL VALUES 
C 

URITE(lOflOl) CALA(I) 

101 FORKATC' '»2E15.5) 

21 CONTINUE 
IF<I0CLOS(10)) STOP 
STOP 

END 
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B.8 LISTING OF INTSTY.FOR 


oRiGfNAL m 
OF POOR QUALITY 

CDCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c c 
c c 
c c 

C PROGRAM TO CALCULATE THE INTENSITY SPECTRUM C 
C C 
C ‘ C 
C C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

c 

C**«**«*«*#m VERSION i 1 ****************** 

C************* PROGRAMMER i R.NAVANEETHAN ****************** 

C««««*#«#***«* DATE : 13-SEP-83 *m***44**m**«« 

c 

C 

C*m444t44#4* A WORD ABOUT THE UAY THIS PROGRAM IS URITTEN !•!!!!!!! 

C444444444444# 

C444444444444# EACH SPECTRAL DATA ARRAY IS ASSOCIATED WITH AN 
C4M4444m##4 ADDITIONAL ARRAY WHICH DEFINES THE RELEVANT PARA 

C4444444444444 METERS ASSOCIATED WITH THE ARRAY, FOR INTENSITY 

C4444444444444 SPECTRUM THESE ARE SAMPLING FREQUENCY (AilALYSIS 

CI444444444444 FREQUENCY RANGE) i MIC SPACING AND AREA ASSOCITAED 
C4444444444444 WITH THE MEASUREMENTS. ADDITIONALLY » ONE MORE SPEC 

C*444444444f44 CAN ALSO BE INCLUDED, THESE FOUR VALUES ARE STORED 

C4444444444444 IN THE BEGINNING OF EACH INTENSITY SPECTRUM DATA 

C4444444444444 ON THE DISK, THE UNIT OF INTENSITY VALUES STORED 

C4444444444444 IS WATT/M-2, 

C 

C234S67 

C 

C444444 DIMENSION STATEMENTS 
C 

COMPLEX AXPS(402)»CAL<402) 

COMPLEX CMPLX 

REAL BBIB(4)»BINT(402)»X(402)»Y(402)fSPL<402) 

REAL SFREQ(100)rAIN<iOO)fBIN<100>fAVG(100)rSPAC(100)fAREAaOO} 

CHARACTER*! ANS 

CHARACTER«15 INAHE(IOO) fNAME 

CHARACTER *3 IN 

CHARACTER *5 OUT 

data 0UT/‘ASP00V 

DATA IN/'AINV 

DATA NAME/'BtXSSOOL,OOOV 

DATA SPL/402*0.0/ 

C 

C#44444 READ MAG CAL FACTORS 
C , 

IF(IOREAD(6f2>0»*A:CALDAT,DAT*)) GOTO 421 
GO TO 422 
421 • CONTINUE 

«RITE(1»423) 

423 FORMAK' ERROR OPENING FILE A:CALDAT,DAT' ) 
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STOP OTG'i\fAl PAGS B 

422 CONTINUE POOR QUAU1T 

READ<6>700) ACALrBCAL 
READ(6r700) ARANfBRAN 
READ(6>701) NAUGAfNAUGB 

700 FQRMAT(2E1545) 

701 F0RMAT(2I5) 

IF (I0CL0S(6)) STOP 
C 

Ct#**** READ XPS FILENAME AND TEST DETAILS 
C 

IF (I0READ<6»2r0»*AtXPSCAT.DAT*)) GOTO 671 
GO TO 672 

671 URITE(1»673> 

673 FORMATS ERROR OPENING FILE AJXPSCAT.DAT') 

STOP 

672 CONTINUE 
J=1 

322 CONTINUE 

READ<6r323>ENDFILE=324)INAHE(J)»SFREQ(J)rAlN(J)fBlN<J)fA06(J) 

SfSPAC(J)fAREA(J) 

J=J+1 
GO TO 322 

323 FORHAT(A12»1X»6FO.O) 

324 JFILE=J-1 
CONTINUE 

IF (I0CL0S(6)) STOP 
C 

C«*t«*4* READ PHASE CAL VALUES 
C 

NAME=*AtCALII.DLO* 

IF (SFREQ(l).GT.lOOl) NAME=*A:CALII,DHI* 

IF<IOREAD(9f2fO»NAME)> GOTO 961 
GOTO 962 

961 WRITE(l»963) 

963 FORMAT(' ERROR OPENING AJCALII.DLO OR AtCALII.DHI FILE') . 

STOP 

962 CONTINUE 

READ<9fl03) AlfA2iA3fA4 
103 F0RHAT(4F15«5) 

DO 120 I=l»402 
READ(9rl21> CAL(I) 

121 F0RMAT(2E15.5) 

120 CONTINUE 

IF (I0CL0S(9>) STOP 
C 

C******* CHANGE DISK FOR OUTPUT FILES 
C 

WRITEdfOOO) 

800 FORMAK' REMOVE PROGRAM DISK IN DRIVE A:» INSERT 
&OUTPUT DISK AND HIT RETURN'!) 

PAUSE 

CONTINUE 


93 


i: 





U i 


ORIGINAL PAGE IS 
IFILE=1 OF POOR QUALITY 

C 

CHANGE N DEPENDING ON THE ANALYSER SPEC 
C 

N=402 

C 

MAIN LOOP FOR FILES 
C 

DO 900 ICsitJFILE 

IF (IFIX(SFREQ(IC)+.05).NE.IFIX<Al+.05)) GO TO 1000 
752 CONTINUE 

IF(lOREAO(9f2»OfINAHE(IC))> 60 TO 7S0 
GO TO 751 

750 CONTINUE 
URITE(lf780) 

780 FORHAT<' CHANGE INPUT FILE DISK IN DRIVE Bt AND HIT RETURN'D 
PAUSE 
60 TO 752 

751 CONTINUE 
WRITE(1^79a) INAHEdO 

798 FORMAT!' 'lAO) 

J=1 

DO 1 I=sl»67 

READ (9»100) (X(Jl)^'Jl=Jf J+5) 

C 

C****** FORMAT STATEMENT DEPENDS ON THE ANALYSER SPEC 
C 

100 FORMAT(6EO.O) 

C 

C#«***0 NEXT STATEMENT DEPENDS ON THE ANALYSER OUTPUT SPEC 
C 

JaJ+6 

1 . CONTINUE 

J»l 

DO 10 I=1t67 

READ(9»100) (Y<Jl)fJlsJ»J+5) 

JsJ+6 

10 CONTINUE 

IF(IOCLOS(9)5 STOP 
BBIB<1)-SFREQ<IC) 

BBIB<2)=2. 

BBIB(3)=SPAC(IC)*»0254 

6BIB<4)=AREA<IC)«*02S4««02S4 

C 

Cfm*# DATA FOR OUTPUT FILE SPEC 
C 

NAHE-INAME(IC) 

CALL PUTCHR(NAMEflfKHAR(IN»l)‘, 

CALL PUTCHR(NAMEf3fKHAR<INf2)) 

CALL PUTCHR(NAME»4fKHAR(INf3)> 

C 

C****** CORRECT FOR HAG AND PHASE CALIBRATION IN DATA CHANNELS 
C 


I 



V i 


; « 


M I 


if ■ 

i . 


\ 


i'l 


it 


m'I, 


•Jj 


S' t 


i* 
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0RIQ(N/\1 i'>j 
OF POOR QUAUTV 

TEHP=ACAL*BCAL*(AIN<IC)*BIM(IC))/(ARAN*BRAN)/<<^M6<IC) 

a/FLOAT(NAMGA)) 

DO 111 I=l»402 
X(I)=XCI)<tTEMP 
Y(I)=Y(I)*TEMP 
AXPS<I)=CMPLX(X(I>»Y<D) 

111 CONTINUE 

DO 4 I=lfN 

AXPS<I)=.AXPS(I)*CAL(I) 

4 CONTINUE 
C 

OPEN a WRITE OUTPUT FILE 
C 

IF(IFILE.LE.41> GOTO /75 
MRITE(1»776) 

776 FORMAK' REMOVE OUTPUT DISK IN DRIVE At AND INSERT 
t INITIALIZED DISK AND HIT RETURN'D 
PAUSE 
IFILE=1 

775 CONTINUE 
758- CONTINUE 

IF(IOWRIT<10»2^0fNAME)) 60 TO 755 
GOTO 756 

755 URITE(1»757) 

757 FORMAT(' CHANGE OUTPUT DISK IN DRIVE A5 AND HIT RETURN'*) 

PAUSE 
GO TO 758 

756 CONTINUE 
IFIL£=IFILE+1 

WRITE(10»123) (BBIB(I)»I=1»4) 

123 F0RMAT(lXf4E15.5) 

C 

C 

C****** CALCULATE INTENSITY 
C 

c 

DO 5 I=i»N 

BINT(I)=AIMAG<AXPS(I)) 

C 

C«**tt* STD SEA LEVEL VALUE FOR DENSITY OF AIR WAS ASSUMED. FOR 
C**##«* BETTER ACCURACY BOTH PRESSURE AND TEMPERATURE MAY BE INPUT 
Ct****4 RHO CALCULATED. 

C 

RH0=1.225 

BINT(I)=BINT(I)/(1.225*BBIB(3)) 

C 

C****** THE FOLLOWING STATEMENTS NEEDS TO CHECKED WITH ANALYSER SPEC 
C 

0MEGA=2.4c3.1415962*FL0AT(I)*BBIBa)*2.56/1024. 

BINT(I)=BINT(I)/OMEGA 

WRITE(10»104) BINT(I) 

104 FORMAT(' '.E15.5) 

5 CONTINUE 
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IF(IOCLOSnO)) STOP OF POOR 

DO 6 I=lf402 

SPLa)=SPL<I)+BINT<I)»BBIB<4) 

6 CONTINUE 
900 CONTINUE 

HAHE=INAHE(1) 

CALL PUTCHR(NAHErl»KHAR(OUTfD) 

CALL PUTCHR(NAHEr3rKHAR(0UT»2)) 

CALL PUTCHR(NAHE»4iKHAR(0UT»3)) 

CALL PUTCHR(NAHE»6fKHAR(0UT>4)) 

CALL PUTCHR(NAME»7»KHAR<0UT»5)) 

MRITE(1»790) 

790 FORHAT(' REHOVE OUTPUT DISK IN DRIVE A*. AMD INSERT SPL DISK'4) 
PAUSE 

CONTINUE 

IF (I0WRIT<9f2»0»NAME)) STOP 
DO 791 I=lf402 

FR=FLOAT (1-1 ) 4SFREQ ( 1 ) *2 . 56/1024 . 

URITE(9f792) FR»SPL(I) 

791 CONTINUE 

792 FORMAT(^ 'f2El5.5) 

IF (I0CL0S(9)) STOP 
GO TO 999 

1000 CONTINUE 

NRITE(lr764) 

764 FORHAK' ERROR IN THE ANALYSIS RANGE SPEC^) 

999 CONTINUE 
STOP 
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ORIGINAL PAGE IS 

C234567 OF POOR QUALir/ 

DIHENSION AINT(402) 

CHARACTERtl2 NAHEfNAHEl 
CHARACTER*! CHRA»CHRB 
HRITE(l»tOO) 

2 CONTINUE 

100 F0RMAT(' ENTER FREQUENCY VALUE OF INTEREST') 

READ(lfiOl) FREQ 

101 F0RMAT(F0»0> 

I CONTINUE 

ICOIJNT = FREQ/1 *25+1 
IF (FREQ*GT.500,05) ICOUNT =FREQ/12.5+1 
IF(FREQ*LE.500.) CHRB = *f 
IF(FREQ.GT. 500,05) CHRB = *H* 

URITE{1,103) 

103 FORHAK' ENTER CATALOG FILE NAME CONTAINING INTENSITY DATA 
JFILE NAMES') 

REA0(lrl04) NAME 

104 FORMAT(AO) 

URITEarl05) 

105 FORMAT<' INSERT FIRST DATA DISK IN DRIVE BJ') 

PAUSE 

3 CONTINUE 
JC=0 

IF (I0READ<6»2»0»NAME)) STOP 

4 CONTINUE 

READ(6»106»ENDFILE=5) NAMEl 

106 FORNAT(AO) 

WRITE(lfl30) NAHEl 

130 FORMAT(' '»A12) 

CALL PUTCHR(CHRArl»KHAR<NAHEl,8)) 

JC-JC+1 

IF(CHRA,NE.CHRB) GO TO 999 
10 CONTINUE 

IF(I0READ<7»2»0»NAME1>) GO TO 8 
GO TO 9 
8 CONTINUE 

URITE<lfl07) 

107 FORHAT(' INSERT 2ND DATA DISK DRIVE IN DRIVE Bt'*) 

PAUSE 

GO TO 10 
9 CONTINUE 

REA0(7fl20) Bl>B2^B3fB4 
120 F0RHAT(4E15,5) 

DO 11 I=1»IC0UNT 
READ<7»108) BINT 

108 F0RMAT(E15.5) 

II . CONTINUE 

IF (BINT.LE.O.) BINT=1,E'-12 
AINT<JC)=10,*(AL0G10(BINT/1,E-12)) 

IF (IQCL0S(7)> STOP 
GO TO 4 

5 CONTINUE 
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IF (I0CL0S(6)) STOP 
JK'sJC 

URITE<1»109> 

109 FORHAir ENTER NAME FOR OUTPUT FILE') 
READdfllO) NAME 

110 FORHAT<AO) 

IF a0URIT<10»2»0»NAHE)) STOP 
DO 12 I-lfJK 
WRITE(lOilll) AINT(I) 

111 FORMAT(' 'fFlStS) 

12 CONTINUE 

IF <IOCLOS(10)) STOP 
GO TO IjOOO 

999 URITE(lf200) 

200 FORMAT(' SOMETHING IS UR0N6 IN FILE NAMESM!') 

1000 CONTINUE 
STOP 
END 
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B.IO LISTING OF ZPLOT.FOR PASS ^3 

OF POOH QUAirn: 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCCCCCCCCCCCCCCC 


c c 

C TRANSHISSION LOSS PLOT PROGRAH C 
C C 
C PROGRAMMER t R. NAVANEETHAN VERSION i ll«i7-83 C 
C C 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C THIS PROGRAH READS THE SOUND POWER DATA FILES TRANSMITTED FROM 

C Z 100 COMPUTER AND PLOTS. THE TL VS FREQUENCY .THIS 

C PROGRAM IS DEVELOPED ON HSNC-li COMPUTER SYSTEM AND HAKES EXTENSIVE 
C USE OF USER WRITTEN (KU-FLIQHT RESEARCH LAB) PLO?'{ER SUBROUTINES. 

C HINC COMPUTER HAS RT-ll OPERATING SYSTEM AND THE FORTRAN USED IS 

C PDP-FORTRAN. MINC HAS EXTENSIVE FORTRAN CALLABLE IEEE-488 

C IB SUB-ROUTINES. THESE SUBROUTINES HAVE BEEN USED IN THIS PROGRAM 

C THE FOLLOWING SUBPROGRAMS ARE REQUIRED*. 

C TLTEST (SOURCE LISTING GIVEN BELOW) 

C OIPLTLIB (WRITTEN AT KU-FRL) 

C OtIBLIB (SUPPLIED BY DEC) 

PROGRAH TLTEST 
C DIMENSION STATEMENTS 

C MAX NUMBER OF POINTS 304 

C TL OR NR = NR(I) 

C FREQUENCY(H2)= FRE<I) 

DIMENSION ARRAY ( 9 ) » ARAY ( 9 ) ? AviRA ( 4 ) » VAX ( 6 ) r FRE < 804 ) 

REAL NR(80«1)jJ>I 
BYTE NAME (20) 

COMMON /ONE/ FRE»NR 
CONTINUE 

DATA ARRAY /'2 '»'4 S'5 'f'6 S'7 'f'8 'f'9 'f 

.'lOOV 

DATA ARAY /'2 S'3 S'4 'f'5 'f'6 't'7 S'8 

,'9 'f'lOOO'/ 

DATA ARRA /'2' f'Z' f'A' t'5' / 

DATA VAX r O'f' 20'»'40'»'60S'80 ’f'100V 
TYPE • 

TYPE If'TURN ON PLOTTERf PLACE PAPER ON PLOTTERi' 

TYPE »r' INSERT DESIRED PEN» RETURN WHEN READY.' 

PAUSE 

CALL PSC(5) 

CALL PCLR 

CALL IBSEND('IP lS43fl488f9559»7520? ' »-l»5) 

CALL SCL (AL0G10(20,)»AL0G10<5000.)f0,fl00.) 

4 TYPE I 

1 FORHAT('<'»'DO YOU WANT AXES DRAWN? (YES=1 N0«2) ') 

ACCEPT *»KL 

CALL CSIZ(1.5»1.5»0»0) 

IF(KL.NE.l) GO TO 21 

CALL PLT (ALOG10(20.)fO.»-2) 

CALL PLT (AL0Gl0(S000.)f0.f2) 

CALL PLT (ALOG10(5000.)flOO.»2) 

CALL PLT (AL0G10(20.)»100.f2) 
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CALL PLT (AL0G10(20.)»0,*-1) 

J » 1 . 

5 DO 10 II=2»iO 

I*IIS£l. 

CALL PLT (AL0G10(U10,*J)»0*f0) 

CALL IPLT (0,f0.f-2) 

IF(I,NE.10) CALL IPLT <0»fl.5r»i) 

IF(I*EQ.10) CALL IPLT (0,»3*r-l) 

IF(I,HE.10)CALL CPLT(-.4»-l .6) 

IF(I,EQ.10)CALL CPLT(-(AL0G10< J4100. )/2. ) r -2.8) 
IF(J.EQ»i) ENCODE (3r6»GN) ARRAY<I>1) 

6 FORMAT ( A3 > 

IF(J.E0»10) ENCODE (4»7»GS) ARAY<I-1) 

7 F0RHAT(A4) 

IF(J.EQ»100) ENCODE (1»8»C0R) ARRA(I-i) 

B FORMAT(Al) 

IF(J,EQ,1) CALL LBL(GM) 

IFU.EQ.IO) CALL LBL(GS» 

IF(J.EQaOO) CALL LBL(COR) 

IF(J.EQ.100 .AND. I.EQ.5) 60 TO 20 
10 CONTINUE 
J * J * 10. 

GO TO 5 

20 CALL PLT (AL0G10(300. ) »0. » 1 > 

CALL IPLT(0.r-6.5fl) 

CALL CPLT(-7.f-1.5) 

CALL CSIZ <2.0fl.5f0»0) 

CALL LBL< 'FREOUENCYf HzM 
CALL PLT <AL0G10(20.>f0.il) 

CALL CSIZ(1.5fl.5r0f0) 

DO 30 IlsOflOOfS 

i=im. 

C PUT TIC MARKS ON VERTICAL AXIS NOU 

CALL PLT(ALOG10(20.)tI*l.rl> 

CALL IPLT<AL0G10(1.07)»0.»2) 

K=7 

IF(I.EQ.O)Kal 
IF(I.EQ.20)K=2 
IF(I.EQ.40>K=3 
IF(I.EQ.60)K=4 
IF(I.EQ.80)K=5 
IF(I.EQ.100)K=4 
IF(K.EQ.7> GO TO 29 
CALL PEN 

CALL aPLT(-4.f-.3) 

VLBL=VAX<K) 

CALL LBL<VLBL) 

29 CALL PLT(ALOG10<20.)»I*l.»l) 

30 CONTINUE 

CALL PLT(AL0G10(20.)»50.,1) 

CALL IPLT<-AL0G10(1.45)f0.fl) 


100 


ORIGINAL PAGE 13 
OF POOR QUALUfY 


CALL CSIZ(2.0»1.5»90.»0) 

CALL CPLK-lO.fOi) 

CALL LBL< 'TRANSMISSION LOSSr dB') 

21 CONTINUE 
TYPE 122 

122 FORMAT!' LON FREQ DATA') 

KIMAXal 
IFLG=»1 

CALL DATRED(KlMAXflFLG) 

KIMAX=KlHAX+0 
IFLGs2 
TYPE 123 

123 FORMAT!' HI FREQ DATA') 

CALL DATRED!KlMAXrIFLG> 

KIHAX=KIMAX-1 

27 CONTINUE 

DO 34 II^i»KIHAX 
‘ X=FRE!II) 

YsNRdl) 

IF!II»EQ,1)CALL PLOTKXfY) 

IF!II*NE.1)CALL PL0T2!XfY) 

34 CONTINUE 
CALL PEN 

TYPE*»'WANT LEAST SQUARE LINE? YES' 

ACCEPT* fILSL 
IF !ILSL .NEd) GO TO 42 
TYPE »» 'ENTER MIN FREQUENCY FOR LEAST SQUARE LINE?' 
ACCEPfltjAHF 

§UMX2=0. 

SUMY2=0* 

N1=0 

DO 118 II=1»KIMAX 

IF!FRE!II) .LT, AMF) GO TO 118 

T1»AL0G10!FRE!ID) 

Nl=NUl 

SUMX=SUHX+T1 

SUMY=SUMY+NR!II) 

SUMX2=SUMX2+T1*»2 

SUHXY=SUMXYm*NR!II) 

118 CONTINUE 

SLaP=!SUMXY-SUMX*SUMY/Nl)/!SUMX2-SUMX*»2/Nl) 

YINT=! SUHY-SLOPKSUMX) /N1 

TYPE*»SLOPfYINT 

X=AL0G10!AHF) 

Y=8L0P*X+YINT 
TYPE*f 'CHANGE PEN ' 

PAUSE 

CALL PLT!X»Y»3) 

X=ALOG10!5000») 
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CALL PLT(XrYr2) 

CALL PEN 
117* nONTlNUE 
42 TYPE 41 

41 rORHAK'fVf'HORE CURVES PLOTTED? <YES*1 N0-2> 
ACCEPT »fHC 
IF<HC.EQ»1) 60 TO 4 
STOP 
END 
C 

SUBROUTINE PLOTKXfY) 

CALL PLT (AL0G10(X)»Y»3) 

RETURN 

END 

C 

SUBROUTINE PL0T2<X»Y> 

CALL PLT (AL0610(X)rYr2) 

RETURN 

END 

C 

SUBROUTINE LABEL 




C 

BYTE LABEL(80) 

DATA LABEL/80* ' V 
C 

TYPE »f 'LABEL DESIRED? YES « 1' 

ACCEPT *rIL 
IF (IL.NE.l) 60 TO 10 

TYPE *» 'ENTER LABEL (LESS THAN 80 CHARACTERS)' 

ACCEPT 5»(LABEL(II)fII«l»80) 

5 F0RNAT<80A1) 

TYPE *f'HOVE PEN TO DESIRED LOCATION OF LABEL' 

TYPE *»'HIT RETURN WHEN SATISFIED' 

PAUSE 

CALL CSIZ (1.5»1.5»0r0) 

CALL LBL( LABEL) 

10 RETURN 
END 
C 
C 

CC 

C 

SUBROUTINE DATRED (KINAXf IFL6) 

COMMON /ONE/ FRE»Tc 
REAL TL<804)»FRE(804) 

BYTE INAMEl(ll)fINAME2(ll) 

87 CONTINUE 
TYPE 222 

222 FORMAK '»'»/»' ENTER THE NAME OF DATA FILE (TEST N/OUT PANEL? 
ACCEPT 23f (INAMEl(II)f Il^lflO) 


') 
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I 


d 


a 




•I 

II 
■ i 

U 


J 


ii 
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ORIGINAL PAQP 

23 FORMAT (10 (AD) OF POOR QUALITY 

OPEN (UNIT=0fNAME=INAHEDTYPE='OLDSACCESS«'SEQUENTIALS 
2 FQRM«'FORMATTED' >ERR=87»REC0RDSIZE=1026) 

307 CONTINUE 
TYPE 322 

322 F0RMAT( '$'»/»' ENTER THE NAME OF DATA FILE (TEST UITH PANEL)? D 
ACCEPT 23f (INAME2(H)fII=l»10) 

OPEN ( UNIT=9 > NAHE=INAME2 » TYPE= 'OLD ' » ACC£SS= ' SEQUENTIAL ' » 

2 F0RH='F0RHATTED'»ERR=387fREC0RDSIZE=1026) 

K=KIMAX 

Kl=sl7* 

K2=401 

IF(IFLG.EG.2) Kl=41 
DO 41 I- D402 
READ(SrlOO) XDYl 
READ(9D00) X2fY2 
100 F0RMAT(2E15,5) 

IF ((I,LT»K1) .0R.(I,GT.K2) ) GO TO 41 
TL(K) = 10!|tALOGlO(Yl)-lO.»ALOGlO(Y2) 

IF(IFLG»EQ.l) FRE(K)= FLOAT( I-l )*1 ,25 
IF(IFLQ.EQ,2) FRE(K)= FL0AT( I-D*12.5 
K=K+1 

41 CONTINUE 

CLOSE (UNIT=8) 

CLOSE (UNIT=9) 

KIHAX=K 

RETURN 

END 


roai 


1 

ll^tl 


»r«K 

I 

^ 't 
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B.ll LISTING OF INTMAP.FOR 


original page b 

OF POOR QUALITY 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 
c 

C INTENSITY MAP INTERPOLATION PROGRAM 

C 
C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C MERSION*. 2 

c programmer: hike graham » r. navaneethan 

C date: 22-NOVEHBER-83 

c 

c 

c 

c 

c dimension statements 

c 

DIMENSION A(9»9)fA,\<9)»AY<9) 

DIMENSION Bl(310)>B2(810)fIB3(810)fIC(810) 

BYTE INAMEa5),0NAMEaS) 

DATA INAMErONAME/30*' • i 
TYPE*» 'ENTER INPUT FILE NAME' 

ACCEPT 122F(INAME(II)fII=lf 14) 

122 F0RMATa4(AD) 

TYPEI» 'ENTER OUTPUT FILE NAME' 

ACCEPT 122»(0NAME<II)fII=1»14) 

123 TYPE»> 'ENTER CONTOUR INTERVAL IN <JB (10. OR 5. OR 2, OR 1.)' 

ACCEPT*, CONINT 

IF( .NOT.((CONINT.EQ.10.>.OR.(CONINT*EQ.5*).OR.(CONINT.EQ.2.) 
l.OR.(CONINT.EQ.l))) GO TO 123 
IF(C0NINT,EQ,10.) IC0N=1 
IF(C0NINT.EQ.5.) ICON = 2 
IF(C0NINT.EQ,2.) .ICON = 3 
IF(CONINT.EQ.l) IC0N=4 

C THE FOLLOWING LOOPS ASSIGN THE VALUES FROM THE INTENSITY TESTS 

C TO AN ARRAY UHICH REPRESENTS THE TEST GRID. 

C 

OPEN (UNIT=8,NAHE=INAME»TYPEs'0LD') 

DO 10 J=lf9 
DO 10 l=lf9 
READ (3i300) A(I>J) 

300 FORMAT (F15.5) 

10 CONTINUE 
C 

C THE FOLLOWING LOOPS TRANSLATE GRID LOCATIONS INTO LOCATIONS 

C DEFINED BY THE DISTANCE FROM THE MPD EDGE. 

C 

DO 11 I=l»9 

AX(I) = 1 + 2.*FL0AT(I-1) 

AYd) = I + 2.*FL0AT(I-1) 
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11 CONTINUE 

CLOSE ( UN IT=8) 

THE FOLLOWING LOOPS DO THE INTERPOLATION BETWEEN POINTS IN 
A HORIZONTAL DIRECTION. THE LOCATION OF EACH DIVISION OF 10 
DECIBELS BETWEEN THE POINTS WILL BE FOUND FOR HAPPING PURPOSES. 

0PEN(UNIT=10fNAHE='INTl .THPSTYPE»'NEH' ) 

DO 20 J=lf9 
DO 20 I=l>8 
Y2 = AUI + l)»J) 

Y1 = A(I»J) 

X2 s AX(I+l) 

XI = AX<I) 

IF(Y2.E0.Y1) GO TO 20 

THESE STATEMENTS ENSURE THAT Y2 IS ALWAYS GREATER THAN Yl. 

IF (Y2.GE.Y1) GO TO 23 
XITEHP = XI 
YITEMP = Yl 
Yl=Y2 
YZaYlTEMP 
X1*X2 
X2=X1TEHP 
23 CONTINUE 

SLOPE = (X2 - X1)/(Y2 - Yl) 

CALL INVAL(Yl,Y2rZCAL.IC0N) 

21 CONTINUE 

IF <ZCAL.6T,Y2) GO TO 22 
XLOC = (SLOPEKZCAL - YD) + XI 
WRITE (lOilOO) XLOC»AY(J),INT(ZCAL) 

ZCAL = ZCAL + CONINT 
GO TO 21 

22 CONTINUE 
20 CONTINUE 

THE FOLLOWING LOOPS DO THE INTERPOLATION BETWEEN POINTS IN A 
VERTICAL DIRECTION. THE LOCATION OF EACH DIVISION OF 10 
DECIBELS BETWEEN POINTS HILL BE FOUND FOR HAPPING PURPOSES. 

DO 30 I=D9 
DO 30 J=lf8 
Y2 = A(D(J+1)> 

Yl = A(I»J) 

X2 = AY(J+1) 

XI = AY(J) 

IF(Y2.E0.Y1) GO TO 30 
C THESE STATEMENTS ENSURE THAT Y2 IS ALWAYS GREATER THAN Yl. 
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IF (Y2.GE.Y1) 60 TO 33 QUALP 

XITEH? - XI 
YITEHP = Y1 
Y1 * Y2 
Y2 s YITEMP 
XI = X2 
X2 = XITEHP 
33 CONTINUE 

SLOPE ■ <X2 - X1)/<Y2 - Yl) 

CALL INVAL<YlrY2»2CAL»ICQN) 

31 CONTINUE 

IF <ZCAL.6T»Y2) GO TO 32 

YLQC = <SL0PE»(ZCA' - YD) + XI 

WRITE (lOflOO) AX(1)»YL0C»INT(ZCAL) 

ZCAL = ZCAL + CONINT 
60 -TO 31 

32 CONTINUE 
30 CONTINUE 

CLOSE{UNIT=10) 

100 FORHAT <lXfF6.2flX»F6.2,lX»I3) 

SORTING OF DATA INTO SEQUENTIAL DIVISIONS OF 10 DECIBELS 

OPEN t UN IT=9 f NAHE* ' INTI » THP ' » TYPE= ' OLD ' ) 

I = 1 

40 CONTINUE 

READ(9»200fEND=45) Bid) rB2( I) »IB3(I) 

200 FORHAT ( lXfF6.2»lXfF6»2f lXrI3) 

I - I + 1 
GO TO 40 

41 CONTINUE 
45 NS = I - 1 

CL0SE(UNIT=9) 

CALL SORT <NSfIB3rIC) 

OPEN ( UN IT= 1 1 » NAHE-ONAHE f TYPE= ' NEW ' ) 

DO 42 J=DNS 
J2 = IC(J) 

URITEdDlOO) Bl(J2)»B2<J2)dB3(J2) 

42 CONTINUE 
STOP 
END 

THE FOLLOWING IS A MODIFIED BUBBLE SORT WRITTEN BY 
R. NAVANEETHAN. 

SUBROUTINE SORT(NSdAf IC) 

DIMENSION IA<810)fKSORT<810)dC(810) 

DO 1 IS=DNS 
IC(IS) = IS 
KSORTdS) s lAdS) 

1 CONTINUE 


Cl W 


ORIGSW-^’*- ^ 
OF POOR QUALITY 


DO 3 

DO 2 JSsl,NS-IS 

IF(KSORT(JS).LE»KSORT(JS+l)> 60 TO 2 
IT = KGORT(JS) 

ITC = IC(JS) 

KSORT(JS) = KSORKJS+l) 

ICUS) = ICUS+1) 

KS0F;T(JS+1> - IT 

IC-IJS+l) = ITC 

CONTINUE 

CONTINUE 

RETURN 

END 


SUSROUTINE INUAL(YlfY2»ZCAL»IC0NT) 

IFdCOHT.NE.l) GO TO 1 
ZCAL=FLOAT(INT(Y1/10»)*10)+10. 

‘30 TO 4 
CONTINUE 

IF(IC0NT.NE.2) GO TO 2 
DEC=(Y1-FLQAT(INT(Y1/10*0)»10))/10, 

IF '.DEC^LT.OtG) ZCAL= FL0AT<1NT(Y1/10, ))*10.+5. 

IF (DEC»GE.0.5) ZCAL = FLOAT! INT( <Yt 4 10 . )/10 . )»10. ) 

60 TO 4 

CONTINUE 

IF(IC0NT.NE.3) 60 TO 5 
IDEC=M0D(INT<Yl)f2) 

IFaifEC.EQ.l) ZCAL=FL0AT(INT<Y1))+1. 

IF(IDEC.EQ.O) ZCAL=FL0AT(INT(Yl)>+2« 

GO TO 4 

CONTINUE 

ZCAL=INT(Y1)+1. 

CONTINUE 

RETURN 

....END 
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CCCCCCCCCCCCCCCCCCCCCCCGCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


t t 

t INTENSITY MAP PROGRAM * 

tt%ttt*tt*t*tt%t*%$ttt******t$tt*t*****t*tnt*ttttttM 


PROGRAMMER : BRIAN QUAYLE 
VERSION 1 11-15-83 

THIS PROGRAM PLOTS INTENSITY MAPS USING DATA GENERATED IN 
THE KU-FRL ACOUSTIC TEST FACILITY ^ 

MORE INFORMATION ABOUT THE PROGRAM CAN BE OBTAINED AT 
THE PROJECT MANAGER OF THE NOISE PROJECT » 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 


* 


PROGRAM PLOT 

BYTE XLBL(39)»INAME(11)»YLBL(37) 

DATA XLBL/'M'r'IS'C'f 'R'f'O'f'P'f'H'r^O'r'NS'E'f^ 'f 


'AS 


'O'l 

IT' 


'R'f 'IS'ZS'OS'NS'TS'AS'LS' S'LS'OS'CS 
'IS 'OS 'NS' 'f f'N'f'C'i'H'f'E' »'S' f')'/ 




DATA YLBL/'M'»'I'f'C'»'R'f'0'»'P'»'H'»'O'»'N'r'E'f' '» 
'V'f'E'f'R'f'T'f 'IS'tS'A'f 'L'f ' '»'L'»'0'»'C'» 
'A'»'T'»'I'»'0'»'N'»' 'r'('f 'I'f 'N'f'C'f 'H'f'E'r'SS')'/ 


TYPE tf' ' 

TYPE *f'TURN ON PLOTTER » LOAD PAPER f AND SET PI AND P2.' 
TYPE *» 'RETURN HHEN READY' 

PAUSE 

CALL PSC(5) 

CALL PCLR 
CALL PEN 


CALL SCL(0.»9.»0,»9.) 

CALL CSIZ(2.il»5»0,i0») 

TYPE*» 'ENTER PEN VELOCITY 1 = 1 CM/SEC r 2 * NORMAL' 
ACCEPTtfIPNVEL 
IFdPNVEL.EQ.DCALL PENVEL 
8 TYPE* f' DO YOU UANT AN AXIS DRAWN? (YES=1)' 

ACCEPT*, AXDRAW 
IF<AXDRAU.NE.1> GOTO 30 
Pl=9, 

P2=9. 

PX3-0. 

PY3=0. 

PX4=9* 

PY4=0. 


f. • :• 
h I 

i I 




i" V 



1 

to j 


n 
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P5=9. OF POOR QUALir**' 

PX6--0 ♦ 

PY6*18, 

PX7*2. 

PY7--=-2» 

P8-1 » 

PX9^1 » 

PY9=-t. 

CALL XAXX(P1»P2»PX3>PX4,P5»PX6»PX7»P8»PX9»XLIIL) 

CALL YAXX(PlfP2fPY3fPY4»P5»PYdfPY7fP8»PY9»YLBL) 

30 IFdPNVEL.ME.DCALL PNVELN 

TYPE »»'WHAT IS THE NAHE OF THE INPUT FILE?' 

ACCEPT 97dINAH‘^(I)fI=i>10> 

97 F0RHAT(10<A1)> 

1267 F0RHAT(lX»F6.2»lXfF6.2»lX»I3> 

OPEN (UNIT=9»NAHE=INAHEfTYPE='0LD' t ACCESS»'SEQUENTIAL ' t 
. FORHs ' FORMATTED ' r ERR:=30 ) 

ISYM=0 

INT=-999 

CALL SCL<0.fl8.rl8,r0.} 

54 READ(9f 1267fEND=57) XPOSfYPOSf INTSTY 
IF<INT,EQ.-999) IVFRST=INTSTY 
IFdNT.NE. INTSTY) ISYH=ISYH+1 
IFdNT.NE. INTSTY) PAUSE 
INT=INTSTY 

, CALL PLT(XP0S>YP0S»1) 

CALL PEN 
CALL SYMB<ISYI1) 

GOTO 54 

57 IVLAST=INTSTY 

TYPE* f 'DO YOU WANT TO LABEL? (YES-1)' 

ACCEPT* r LAB 

IF(LAB.EQ.l) CALL IDENTCdVFRSTf IVLASTrlSYH) 

467 TYPEIf'DO YOU WANT TO PLOT ANOTHER HAP? (YES*l)' 

ACCEPT* » MAP 
IF(MAP.EQ.l) GOTO 8 
STOP 
END 
C 
C 
C 

SUBROUTINE 1DENTC(IVFRST»IVLAST»ISYH) 

DIMENSION NUM(3) 

BYTE NAME(81)fNUM 
DATA NAME/81*' '/ 

DATA NUM/3*' '/ 

100 TYPE *» 'ENTER LABEL (LESS THAN 80 CHARACTERS)' 

TYPE *» 'ENTER 'ZZZ* IF FINISHED' 

ACCEPT 30f(NAME(J)f J=l»80) 

IF((NAME(1)«EQ»'Z').AND,(NAME(2).EU.'Z').AND.(NAME(3)«E0.'Z')) 
.GOTO 44 
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TYPE *f' ' 

TYPE *f'USE PLOTTER PEN CONTROL TO POSITION LABEL' 
TYPE *f'RETURN WHEN SATISFIED' 

PAUSE 

CALL CSIZ(1.5fli5»O»f0> 

20 CALL LBL(NAHE) 

DO 50 I=l»80 
J=81-I 

IF<NAME(J) .NE.' ') GOTO 60 
50 CONTINUE 
60 CONTINUE 

CALL CPLT(anAT<-J)f-1.2) 

GOTO too 

66 INCRHT=aVLAST-IVFRST)/(ISYM-l) 
IUFRST=IVFRST-INCRHT 
CALL CPLT(3f-2.2) 

DO 112 I=lfISYM 
IVal 

INUM=IMFRST+INCRHT 

IVAR=INUH 

ENC0DE(3fl3fNUH) INUH 
13 F0RHAT(I3) 

CALL SYMB(IV) 

CALL CPLT(1».4) 

CALL LBL('=5') 

CALL CPLTdrO) 

CALL LBL(NUH) 

CALL CPLTdfO) 

CALL LBL('DB') 

IF(IVAR»LT.100) CALL CPLT(-9»-l .2) 

IF<IVAR*GEdOO) CALL CPLT(-10»-l .2) 

112 CONTINUE 

CALL CSIZ(2.5d*5»0.f0) 

30 FORMAT (80A1) 

RETURN 

END 

C 

C 

C 

SUBROUTINE XAXX (PI »P2»P3»P4»P5»P6»P7»P8f P9fLABEL) 

COMMON /ADORES/ INSTR 

COMMON /XYSCLP/ SCLPlfSCLP2iSCLP3fSCLP4 

BYTE ARRAY (10 )» LABEL (3?) 

DATA ARRAY /lO*' '/ 

P15--.1*P9 

P11=»1*P9 

P16=0. 

SCLPi=P6-P7*P3 

SCLP2=P6+(P1“P3)*P7 

HGHT=40./P2 

CALL CSIZ(2.il.5fO,fO.) 
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CALL 3CL<0f »P1»0.»P2) OF POOR QUALITY 

CALL PLT(Plf0.r-2) 

CALL PLT(P1»P2»-U 

CALL PEN 

XPLT=P3 

YPLT==P4-P11 

CALL PLT<XPLTfYPLT»l) 

JINT=0 

IB3=0 

CALL FINT6R(P6tP7»JINT) 

IF(JINT.EQ,l)GOTQ 201 

CALL FINTGR(P6»10. »P7*10f » JINT) 

JINT=JINT+i 

IG3=1 

201 CONTINUE 

IF( (P6+P7*P16) .GT. 10000. ) JINT=0 
4 CONTINUE 

IF(JINT.EQ.O) GOTO 202 
15 F0RMAT<1I9) 

IF(JINT.EQ.1)ENC0DE<9»15»ARRAY(D) IMT(P6+P16*P7) 
IF(JIHT.EQ.2)ENC0DE(9fl.2rARRAY(l)) P6+P16*P7 
G0T043 

202 CONTINUE 

10 F0RHAT<1PE9,2) 

11 FORMAT(1F9.0) 

12 F0RMAT(1F9.1) 

13 F0RHAT(1F9.2) 

14 F0RHAT(1F9.3) 

ENCODE (9»10»ARRAY(1>) P6+P16IP7 
IF(ABS(SCLP1),LT,10000.AND.ABS<SCLP2).LT-^, 10000) 
.ENCODE (9.UfARRAY(l)) P6+P16»P7 
IF(ABS(SCLP1).LT. 1000. AND. ABS<SCLR2) .LT.IOOO) 

.ENCODE (9»12rARRAY(D) P6+P16IP7 
IF(ABS(SCLP1) .LT.100,AND.ABS(SCLP2).LT.100) 

.ENCODE (9»13»ARRAY(D) P6+P16IP7 
IF(ABS(SCLP1) .LT,10.AND.ABS<SCLP2).LT.10) 

.ENCODE (9»14»ARRAY(D) P6+P16»P7 

43 CONTINUE 

IF(ARRAY(1) .NE.' ') GOTO 44 
DO 45 IC=1»9 
ARRAY(IC)=ARRAYaC+l) 

45 CONTINUE 

ARRAYaO)= ' ' 

GOTO 43 

44 CONTINUE 
ISA0=9 

DO 666 IA=1»9 

IF(ARRAY(10->IA).EQ.' ') GOTO 666 
XSAU=10-IA 
UDTO 667 
666 CONTINUE 
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667 CONTINUE 
P1S*ISAM 
DO 888 JJ=1»9 

•• IF(ARRAY(JJ).EQ.'0O ARRAY' 

888 CONTINUE 
P17=0. 

CALL PEN 

CU=-<P18/2.+»25) 

CH=*75*P9+»2 
CALL CPLT(CW»CH) 

CALL LBL( ARRAY) 

C«*-<P18/2.-»25) 

CH*-.2-*75*P9 
CALL CPLT(CMfCH) 

3 CONTINUE 
P16SP16+1. 

P17spi7+1. 

IF<P16tLE»P5) CALL IPLT(0 » fPII »2) 

IF<P16,LE.P5) CALL IPLTd . »0, r2) 

IFCP16.LE.P5) CALL IPLT(0 * fP15»2) 

IF(P17aT.P8) GOTO 3 
IF(P16.LE.P5) GOTO 4 
DELX--*5«P5 
DELY=0. 

CALL IPLT(DELX»DELYfl) 

DO 5 I=lf39 

IF(LABEL(40-I).E0.' ') GOTO 5 
ISAVE=40-I 
GOTO 6 
CONTINUE 
CONTINUE 

P17*FL0AT(ISAVE> 

CALL CSIZ(2.5»l,5f0.»0,) 

CUS-P17/2. 

CH=(1.75*P9-,25) 

CALL CPLT(CUfCH) 

CALL LBL< LABEL) 

CALL PEN 

I F ( SCLPl . LT . SCLP2 . AND . SCLP3 . LT , SCLP4 ) 

» C ALL SCL ( SCLP 1 f 3CLP2 r SCLP3 f SCLP4 ) 

RETURN 
END 
C 
C 
C 

SUBROUTINE YAXXCPl fP2»P3»P4»P5»P6»P7»P8»P9f YLABEL) 

COMMON /ADDRES/ INSTR 

COMMON /XYSCLP/ SCLPl »SCLP2»SCLP3f SCLP4 

BYTE ARRAY(10)»YLABEL(37) 

DATA ARRAY /lOt' V 
P10=.UP9 
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P14*-*1»P9 
PI 9*0. 

SCLP3*P6-P7»P4 
SCLP4*P6+(P2-P4)*P7 
CALL CSIZ(2.il»5^0.f0.> 

CALL SCL(0.rPlr0,»P2) 

CALL PLT(PlfO.»-2) 

CALL PLT(O.fO.r-l) 

CALL PEH 

XPLT*P3-P10 

YPtT=P4 

CALL PLT(XPLT>YPLTfl) 

JINT=0 

IG3=0 

CALL rrHTGR(P6»P7f JINT) 

IF(JINT.EQ.1)G0T0 201 

CALL FINTGR(P6*10.»P74(10.f JINT) 

JIHT*JINT+1 

163=1 

CONTINUE 

IF<<P6+P7*P16),GT, 10000. >JINT=0 
4 CONTINUE 

IF<JINT.EQ,0) 6070 202 
F0RHAT(1I9) 

IF(JINT.EQ,l)ENC0DE(9fl5»ARRAY<l)> INT(P6+P14*P7) 

IF ( JINT . EQ . 2 ) ENCODE ( 9 r 12 f ARRAY ( I > ) P6+P16*P7 

G0TQ43 

CONTINUE 

10 F0RMAT(1PE9.2) 

11 F0RHAT(1F9.0) 

12 F0RHAT(1F9.1) 

13 F0RMAT(IF9.2) 

14 F0RHAT(1F9.3) 

ENCODE (9flO»ARRAY<D) P6+PU»P7 
IF ( ABS(SCLP3). UT. 10000. AND. ABS(SCLP4).LT. 10000) 
.ENCODE (9fll?fiRRAY<l>) P6+PU*P7 
IF(ABS(SCLP3).LT.1000.AND.ABS<SCLP4).LT.1000) 

.ENCODE (9»12»ARRAY(1)) P6+P16»P7 
IF(ABS(SCLP3).LT.100.AND.ABS(SCLP4).LT.100) 

.ENCODE (9fl3»ARRAY(D) P6+P16*P7 
IF(ABS(SCLP3).LT.10.AND.ABS<SCLP4).LT.10) 

.ENCODE (9»14»ARRAYa>) P6+P16»P7 

43 CONTINUE 

IF<ARRAY<1).NE.' ') GOTO 44 
DO 45 IC=1»9 
ARRAY(IC)=ARRAY(IC+1) 

45 CONTINUE 

ARRAY (10)*" " 

GOTO 43 

44 CONTINUE 
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ISAVs? 

DO 666 IA»lf9 

IF(ARRAY(10-IA).EQ.^ ') GOTO 666 
' ISAV-IO-IA 
GOTO 667 

666 CONTINUE 

667 CONTINUE 
P18*ISAV 

IF(P13,GT,P19) P19=P18 
DO 888 JJ-lf9 

IF(ARRAY<JJ).E0.'0') ARRAY( JJ)»'0' 

883 CONTINUE 
P17*0. 

CALL PEN 

CUs<P9f<P18/2,+2.)-Pl8/2.) 

CH=-,3 

CALL CPLT<CU»CH) 

CALL LBL(ARRAY) 

CU»~<P18/2*+P9*<P18/2.+2.)) 

CH=.3 

CALL CPLT(C«fCH) 

3 CONTINUE 
P16=P16+1. 

P17=P17+1. 

IF(P16.LE»P5) CALL IPLT<P10»0, »2) 

IF<P16.LE.P5) CALL IPLT(0, f 1. i2) 

IF(P16.LE.P5) CALL IPLT(P14»0. »2) 

IF(P17.LT.P8) GOTO 3 
IF(P16.LE.P5) GOTO 4 
DELX^iO. 

DELY--.54P5 

CALL IPLT(DELXrDELY»l) 

DO 5 l3l,37 

IF<YLABEL(38-I).EQ.' ') GOTO 5 

ISAyE=38-I 

GOTO 6 

5 CONTINUE 

6 CONTINUE 
P17-FL0AT(ISAME) 

CW=P9*<P19+2.)+»5 

CRsO . 

CALL CPLT(CW»CH) 

CALL CPLT(-1«»0*) 

CALL CSIZ(2,5»1.5.»90.»0.) 

CU»-P17/2. 

CH=0* 

CALL CPLTtCU»CH) 

CALL LBLrrLABEL) 

I F ( SCLPl » LT . SCLP2 , AND . SCLP3 . LT , SCLP4 ) 

.CALL SCL(SGLPlfSCLP2fSCLP3»SCLP4> 

RETURN 
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END 


SUBROUTINE SYMB(ISYH) 

COMHON/ADDRES/INSTR 
BYTE HESSAG(80) 

DATA «ESSAG/80*0/ 

CALL CSIZ(2.5fl.5»0.f0.) 
ENC0DE(l2f2421fHESSAG(i))'SI.175,.35> ' 

2421 F0RHAT(1A12) 

CALL IBSEND<HESSAGE,11»INSTR) 

IFdSYH.NE.DGOTO 1 

ENC0DE(42»2424»MESSAG<l))'UC-99»-3»-3f99»6»0»0f6r-6fOrOf-6» 
l-99f3f3J ' 

2424 F0RHAT(1A42) 

CALL IBSEND<HESSAGf41»INSTR> 

GOTO 99 

1 CONTINUE 
IF(ISY}1*NE.2)60T02 

ENCODE(38f2425»HESSAG<l))'UC-99rOr4f99»-3f-6»6fO»-3»6»-99» 
I0r-4f ' 

2425 F0RHAT(1A38) 

CALL IBSENB(HESSAG»37»INSTR) 

GO TO 99 

2 CONTINUE 
IF(ISYH.NE.3)60T03 

ENC0DE(39»2426rHESSAG<l))'UC-99»-3»2»99rA»0»-3»-6»-3»6»“99» 
l3f-2» ' 

2426 F0RHAT(1A39) 

CALL IBSEND<HESSAG»38fINSTR) 

GO TO 99 

3 CONTINUE 
IF(ISYM.NE.4)60T0 4 

ENC0DE(62f2427»HESSAG(l))'UC-99»-lf3 ''9r-2»-2»0»-2»2»-2f 2f 
I0r2»2»0»2f-2f2»“2f0r-99fl»-3f ' 

2427 FORHATdA62) 

CALL IBSEND<HESSAGf61»INSTR) 

GOTO 99 

4 CONTINUE 
IF<IsyH.NE.5)G0T0 5 

ENC0DE(43»2428jMESSAG(l))'UC-99»3»0»99,-3f4f-3»-4»3»-4f3f 
tAf-99r-2fOf ' , 

2428 F0RHAT(1A43) 

CALL IBSEND<MESSAG»42dNSTR) 

GOTO 99 

5 CONTINUE 
IF(ISYM.NE.6)60T0 A 

ENC0DE(38F2429fHESSA6(l) >'UC-99f2»2»99f-6»0»6»-6»0»6» 
J-99r-2f-2J ' 

2429 F0RMAT(1A38) 
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CALL IBSEND(HESSA6r37fINSTR> 
GOTO 99 
6 CONTINUE 

99 CONTINUE 

RETURN 
END 
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SUBROUTINE PENVEL 

COHNON/ADDRES/INSTR 

CALL IBSENO<^MSir ^-IrlNSTR) 

RETURN 

END 


SUBROUTINE PNVELN 

COHHON/ADDRES/INSTR 

CALL IBSEND<'VS;'»-i»INSTR) 

RETURN 

END 


SUBROUTINE FINT6R(G1 f02» JINT) 

JINT=0 

IF< ( (61-INT<6D) »EQ*0) .AND. < (62>INT(62) ) .EO.O)) JINT»1 

RETURN 

END 
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